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SEMICONDUCTOR  DEVICE  CONCEPTS 
Abstract 

The  effects  of  uniaxial  compression  and  of  hydrostatic  pressure 

on  the  direct  and  indirect  tunneling  processes  in  germanium  tunnel  diodes 

have  been  studied  experimentally  under  forward  and  reverse  bias  at  4.2*K 

and  compared  with  Kane's  theory.  The  diodes  were  formed  by  alloying 

indium  doped  with  3/8  weight  per  cent  gallium  on  (100)  and  (110)  faces 

18  3 

of  germanium  bars  containing  an  antimony  concentration  of  5.5  x  10  /cm  . 
The  first  order  change  of  the  tunneling  current  with  stress  was  measured 
at  fixed  bias  voltages.  For  biases  smaller  than  8  rav  the  current  is 
direct  and  not  affected  by  the  relative  shifts  of  the  (111)  conduction 
band  valleys.  In  the  bias  range  of  Indirect  tunneling  the  anisotropic 
tunneling  from  the  (111)  valleys  was  observed  in  agreement  with  theory. 

In  the  range  of  direct  tunneling  to  the  (000)  conduction  band,  the 
current  change  is  correlated  with  the  stress  Induced  change  of  the  direct 
band  gap  and  of  the  energy  separation  between  the  (111)  and  (000) 
conduction  bands.  This  separation  was  found  to  be  0.160  ±  0.005  ev  ft 
sero  stress  in  agreement  with  optical  measurements  on  degenerate 
germanium.  Some  details  of  the  bias  dependence  of  the  pressure  effect 
including  some  fine  structure  at  small  biases  remain  unexplained. 

Electrical  transport  measurement  have  been  made  on  p-type  ZnTe  and 
n-type  ZnSe.  In  ZnTe  crystals  doping  with  Cu,  Ag  and  Au  produces 
acceptor  levels  at  0.15,  0.11,  and  0.22  ev  respectively.  An  acceptor 
with  an  ionisation  energy  of  0.048  ev  was  found  in  the  undoped  crystals 
and  is  identified  as  the  first  charge  state  of  the  Zn- vacancy.  A  shallow 
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donor  state,  at  approximately  0.01  ev  below  the  conduction  band, 
was  found  in  n-type  ZnSe.  It  also  proved  possible  to  prepare 
degenerate  ZnSe.  The  scattering  mechanisms  limiting  the  lattice 
mobilities  of  both  materials  were  considered.  It  was  found  that 
the  polar  interaction  with  the  longitudinal  optical  phonons 
dominates  the  scattering  of  electrons  in  ZnSe.  This  mechanism 
probably  also  predominates  in  the  scattering  of  the  holes  in  ZnTe. 
However,  the  non-polar  interaction  with  the  optical  modes  could 
also  contribute  significantly  if  the  appropriate  coupling  parameter 
Is  larger  than  we  presently  believe. 

Luminescence  from  p-n  junctions  constructed  from  the  mixed 
crystal  GaAs.  P  and  from  GaSb  has  been  studied.  Visible  red 
radiation  is  produced  by  the  former,  but  at  present  the  efficiency 
is  poor  and  the  junction  profile  is  quite  irregular,  presumably 
because  of  poor  crystal  structure.  Attempts  to  produce  coherent 
radiation  from  these  junctions  has  thus  far  been  unsuccessful. 
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SEMICONDUCTOR  DE/ICE  CONCEPTS 

A.  The  Effect  of  Elastic  Strain  on  Interband  Tunneling  in  Sb-Doped 
Germanium  (H.  Fritzsche  and  J.  J.  Tlemann) 

The  research  reported  here  evolved  from  experiments  carried  out 
under  Contract  AF  19(604) -6623  which  preceded  the  present  one.  These 
early  investigations  were  described  in  Reports  3A  and  4A  of  that  contract, 
dated  23  December  1960  and  28  March  1961.  In  view  of  this  history  it 
seems  appropriate  to  present  here  a  report  of  this  work  which  has  now 
reached  a  satisfactory  conclusion. 

I.  INTRODUCTION 

As  a  result  of  extensive  experimental  and  theoretical  efforts, 

the  main  features  of  the  tunneling  process  in  Esaki  tunnel  diodes^* 

have  been  clearly  established.  However,  since  it  is  difficult  to  assess 

the  validity  of  some  of  the  simplifying  assumptions  and  approximations 

which  underlie  our  present  theoretical  understanding  of  this  process, 

2  3 

it  is  not  clear  to  what  extent  theexisting  theories  *  should  explain 
the  finer  details  of  the  experimental  observations.  This  probelm  is 
particularly  difficult  to  resolve  because  tunnel  diodes  can  only  be 
made  in  highly  impure  materials.  This  fact  precludes  the  possibility  of 
ever  performing  a  tunnel  diode  experiment  under  very  ideal  conditional. 

It  is  hoped  that  some  of  these  difficulties  can  be  alleviated  by 
performing  a  series  of  experiments  upon  each  diode.  In  this  paper  we 
discuss  measurements  of  the  effect  of  hydrostatic  pressure  and  of 
uniaxial  compression  on  the  tunneling  current  of  germanium  tunnel  diodes 
at  helium  temperatures.  A  subsequent  peper  will  report  the  temperature 
dependence  of  the  tunneling  current  of  the  seme  diodes. 
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The  effects  of  elastic  strain  will  strongly  depend  on  the 
effective  mass  anisotropies  and  on  the  deformation  potentials  of  the 

4 

relevant  band  extrema.  Most  of  these  quantities  are  known  for  germanium. 

Furthermore,  in  germanium,  two  types  of  tunneling  have  been  observed,^’** 

3 

and  theoretical  expressions  for  both  types  have  been  derived.  These 
two  tunneling  processes  are  (1)  direct  tunneling  and  (2)  indirect  or 
phonon  assisted  tunneling.  The  relative  proportion  of  the  current 
carried  by  these  two  processes  depends  on  the  temperature,  the  bias 
voltage,  and  on  whether  As,  P  or  Sb  is  used  as  the  donor  impurity.^ 
Although  the  mechanism  is  not  understood  theoretically  at  present, 
direct  tunneling  predominates  in  As-  and  P -doped  diodes  at  all  temperatures 
and  biases.  In  Sb-doped  diodes,  on  the  other  hand,  one  observes  at 
helium  temperatures  several  distinct  bias  regions  in  which  one  or  the 
other  of  these  processes  is  dominant.  This  material  has  the  advantage 
that  both  direct  and  indirect  tunneling  are  observable  in  separate 
bias  ranges  in  a  single  diode. 

This  is  important  because  there  is  considerable  experimental 
uncertainty  as  regards  the  details  of  the  impurity  distribution  within 
the  junction  Itself.  By  measuring  both  types  of  tunneling  in  one 
sample,  this  uncertainty  does  not  affect  the  conclusions.  Antimony- doped 
germanium  tunnel  diodes  have  the  additional  advantage  that  the  unambiguous¬ 
ly  direct  tunneling  to  the  (000)  minimum  in  the  conduction  band  is  much 
more  clearly  observed. 

The  main  effect  of  elastic  strain  on  the  tunneling  current  arises 
from  the  shifts  in  energy  of  the  conduction  band  valleys  with  respect 
to  the  valence  band.  Hydrostatic  pressure  causes  all  of  the  conduction 


band  valleys  to  rise  in  energy,  but  the  (000)  valley  moves  faster  than 

a 

the  (111)  valleys.  In  addition,  shear  stress  can  shift  the  (111) 

o 

valleys  with  respect  to  one  another.  These  effects  can  be  quantitatively 

calculated  by  applying  deformation  potential  theory  to  the  expressions 

3 

for  the  tunnel  current  derived  by  Kane.  The  calculated  behavior  is  in 

good  qualitative  agreement  with  the  main  features  of  the  experimental  data. 

II.  THEORETICAL  CONSIDERATIONS 
3 

According  to  Kane's  theory,  the  tunneling  current  flowing  from  a 
single  valley  in  the  conduction  band  to  the  valence  band  is  given  by 

Id  -  Cd  D  exp  (-00  with  «-  ^dEg3/2  mn^2/F  (1) 

when  the  transition  is  direct,  i.e.  when  the  conduction  and  valence  band 
extreoia  occur  at  the  same  wave  number  k.  It  is  given  by 

Ii  -  Ct  D  exp  (-f)  with  p  -  Xt  E  3/2m*1/2/F  (2) 

when  the  transition  is  indirect,  i.e.  when  the  band  extreme  occur  at 
different  k. 

In  Eqs.  (1)  and  (2)  E  is  the  appropriate  band  gap,  m*  - 

o 

(l/m^  +  1/m^)"1  is  the  reduced  effective  mass  in  the  tunneling  direction. 
The  quantities  C  and  X  are  given  by  Kane.  The  average  junction  field  is 
given  by 

F  -  (2 TTn*  (Eg  +  £n  +  f p  -eV)/K]1/2  (3) 

h 

where  n  -  np/ (n+p)  is  the  reduced  doping  constant,  ft  is  the  dielectric 
constant,  £  n  and  are  the  Fermi  level  penetrations  into  the  conduction 
band  and  the  valence  band,  respectively,  and  E  is  the  smallest  band  gap. 

In  Eq.  (2)  the  energy  of  the  phonon  which  is  required  to  conserve  wave 
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number  in  the  indirect  transition  has  been  neglected  in  the  exponential 
factor.  The  phonon  energy  has  to  be  considered  in  the  density  of  states 
factor  D,  however,  as  will  be  discussed  later.  The  density  of  states 
factor  is  given  in  Kane's  notation  by 

D  -  J[1  -  exp  <-2Es/Ea)]  [fl  (E1>  -  f2  (E2)]  dE.  (4) 

When  there  is  more  than  one  valley  in  the  conduction  band,  or  more 
than  one  valence  band,  the  total  current  will  be  a  sum  over  all  combinations 
of  conduction  band  and  valence  band  extrema  which  can  contribute.  Because 
of  the  exponential  dependence  on  the  reduced  mass  appearing  in  Eqs.  (1) 
and  (2),  however,  this  sum  for  germanium  can  be  reduced  to  only  those 
transitions  which  involve  the  light  hole  band.  Furthermore,  since  the 
conduction  band  valleys  in  germanium  are  quite  anisotropic,  and 
since  the  tunneling  probability  depends  only  on  the  effective  mass 
along  the  direction  of  tunneling1**  the  current  contributions  of  the 
four  (111)  valleys  will  in  general  be  different. 

That  part  of  the  applied  stress  which  is  associated  with  a  volume 

change  affects  the  tunnel  currents  primarily  through  E  and  m*,  which 

© 

appear  in  the  exponents  of  Eqs.  (1)  and  (2).  The  variations  of  these 
quantities  with  stress  can  be  estimated  from  the  known  pressure 
coefficients  of  the  band  gaps. 

The  pure  shear  part  of  the  stress  causes  the  conduction  band 
valleys  to  shift  in  energy  with  respect  to  one  another  in  such  a  way 
that  their  average  energy  is  unchanged.  Therefore,  for  the  shifts  in 
valley  energies  to  produce  a  current  change  which  is  linear  in  stress, 
it  is  necessary  (a)  that  the  degeneracy  of  the  valleys  be  removed  by 
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■hear,  and  (b)  that  the  current  contributions  of  the  valleys  which  are 
raised  in  energy  be  unequal  to  the  contributions  from  those  which  are 
lowered. 

If  the  electric  field  is  in  a  [100]  direction,  all  (111)  valleys 
have  the  same  effective  mass  component  in  the  field  direction.  For 
such  a  diode  an  energy  shift  of  the  valleys  should  not  produce  a  current 
change  which  is  linear  in  shear  stress. If  the  electric  field  is  in 
the  [110]  direction,  however,  two  valleys  (those  along  [111]  and  [ill]) 

have  reduced  effective  masses  along  the  field  direction  m.  ■  m_  ■ 

-1  1  1 
1  i  L 

(  — —  +  — )  ■  0.027  m  ;  while  the  two  other  valleys  (along  [111]  and 

m  ° 

n  ^  ^  j  2  11 

[111])  have  reduced  effective  masses  m,  -  m.  ■  (r - +  ■= - +  — )  ’  m 

3  4  3mte  3m/e 

0.034  mQ,  where  mte,  m ,  ,  and  are  the  principal  electron  masses  and 
the  light  hole  mass,  respectively.  Even  though  these  effective  masses 
differ  by  only  about  25%,  the  currents  themselves  will  differ  by  a 
large  factor  because  of  the  exponential  dependence.  In  our  samples, 
for  example,  the  exponent  is  about  equal  to  16.  This  leads  to  a 
factor  of  6  difference  in  the  currents.  Such  a  diode  will  exhibit 
a  first  order  stress  coefficient  for  a  shear  which  causes  the  similar 
valleys  to  move  in  the  same  direction.  (On  the  other  hand,  a  shear  which 
raises  one  valley  of  each  pair  while  lowering  the  other  valley  will  have 
no  first  order  shear  dependence.) 

III.  EXPERIMENTAL 

Our  experimental  diodes  were  chosen  so  that  they  are  as  identical 
as  possible  except  that  one  sample  has  the  electric  field  along  the 
[001]  axis  and  the  other  along  [110].  For  both  samples  the  uniaxial 
■tress  direction  was  [110].  The  samples  were  in  the  form  of  rectangular 
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bars  with  the  long  dimension  along  [llo]  and  the  short  dimensions  along 
[110]  and  [001].  Two  diodes  were  alloyed  near  the  center  of  each  bar 
on  opposite  faces.  It  was  found  to  be  necessary  to  do  this  because  it 
was  impossible  to  avoid  flexing  the  bars  slightly.  Since  flexure  causes 
equal  and  opposite  stresses  on  opposite  faces,  the  average  of  the  results 
obtained  from  these  diodes  is  independent  of  flexure. 

The  starting  material  was  pulled  along  the  [110]  axis  from  a 
melt  containing  6%  Sb  by  weight.  After  orientation  by  x-rays,  the 
crystal  was  cut  into  wafers.  These  wafers  were  reoriented  and  cut 
into  bars  with  the  long  dimension  along  [llo]  and  the  short  dimensions 
along  [110]  and  [001].  From  resistivity  and  Hall  effect  measurements, 
the  Fermi  level  penetration  of  our  samples  was  calculated  to  be 
.020  ±  .002  volts.  The  diodes  were  formed  by  alloying  Indium  dots  doped 
with  3/8  weight  X  gallium  on  the  appropriate  faces.  Double  contact 
wires  were  attached  to  the  dots  and  the  samples  were  etched  to  remove  the 
perimeter  of  the  junctions.  (This  was  done  because  the  perimeter  has 
the  wrong  orientation.)  The  diameter  of  a  typical  diode  dot  was  about 
0.05  cm. 

The  Fermi  level  on  the  p-type  side  was  estimated  to  be  approximately 
.140  ±  .020  ev.  No  direct  measurements  were  made  to  confirm  this  estimate, 
but  since  none  of  our  conclusions  is  sensitive  to  the  choice  of  as  long 
as  >  £n,  the  accuracy  of  the  estimate  is  unimportant. 

Uniaxial  compresslonal  stresses  X  varying  between  5  x  107  and 
5  x  10  dynes/cm  were  applied  parallel  to  the  [110]  direction  at 
temperatures  between  1.5  and  4.2*K.  The  stress  tunneling  coefficient 
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defined  as  ir  -AI/IX  averaged  over  the  two  opposed  diodes  was  found  to 

be  Independent  of  stress  and  of  temperature  within  these  ranges.  The 

stress  tunneling  coefficients  were  measured  at  fixed  bias  voltages. 

In  order  to  minimize  the  effect  of  series  resistances  due  to  the  leads 

and  the  sample,  the  bias  was  measured  potentlometrlcally  using  one  pair 

of  the  double  contact  wires  attached  to  the  diode  and  to  an  ohmic  contact 

In  close  proximity  of  the  diodes  while  the  current  was  passed  through 

the  other  pair  of  contact  wires.  Liquid  helium  was  used  as  the  pressure 

transmitting  fluid  for  the  hydrostatic  pressure  measurements.  These 

3 

were  extended  up  to  p  ■  10  psi. 

Fig.  1  shows  the  stress  tunneling  coefficient  as  a  function  of 
bias  for  the  two  samples.  In  sample  1  the  diodes  are  on  the  (001)  faces. 
This  sample  should,  therefore,  not  respond  to  the  shear  Induced  shifts 
of  the  conduction  band  valleys.  In  sample  2  the  diodes  are  on  the  (110) 
faces.  The  difference  between  these  two  curves  Is  due  mainly  to  the 
nonequivalence  of  the  two  pairs  of  valleys  for  tunneling  In  the  [110] 
direction  and  the  fact  that  the  shear  part  of  the  applied  stress  causes 
these  pairs  to  move  In  opposite  directions. 

The  same  figure  shows  the  hydrostatic  pressure  coefficient  of 
sample  1  at  4.2*K.  In  order  to  make  a  direct  comparison  with  the 
uniaxial  compression  data  possible  ir^«AI/3Ip  was  plotted.  The 
difference  between  the  uniaxial  and  the  hydrostatic  pressure  coefficients 
of  sample  1  Is  due  to  the  shear -Induced  effective  mass  changes  of  the 
electrons  and  light  holes  In  the  tunneling  direction. 
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The  I-V  characteristic  of  one  of  the  diodes  of  sample  2  is 
plotted  on  Fig.  2.  The  two  diodes  were  so  closely  matched  that 
their  characteristics  did  not  differ  by  more  than  l0°/o  over  the 

whole  voltage  range.  The  characteristic  has  the  shape  typical  for 

5 

indirect  tunneling.  Only  a  very  small  current  can  flow  until  the 
bias  is  large  enough  to  permit  the  emission  of  a  low  energy  phonon 
needed  for  wave  number  conservation  in  the  tunneling  process  from 
a  (111)  conduction  band  valley  to  the  (000)  valence  band.  The 
current  increases  again  when  a  higher  energy  phonon  of  the  same 
wave  number  can  be  emitted.  The  phonons  have  been  identified 
previously^  as  the  TA  and  LA  [111]  phonons  which  have  the  energies 
0.0076  eV  and  0.028  eV,  respectively.*-2 

The  I-V  characteristic  of  the  same  diode  is  plotted  on  a 
different  scale  in  Fig.  3  in  order  to  show  the  rapid  Increase 
of  the  reverse  current  at  a  bias  voltage  of  about  -140  nW.  This 
so-called  Kane  kink  occurs^  when  the  back  bias  is  large  enough 
to  move  the  higher  lying  conduction  band  edge  at  k  ■  0  on  the 
n-type  side  below  the  Fermi  level  in  the  valence  band  on  the 
p-type  side  and  thus  to  allow  a  large  direct  tunneling  current 
to  flow.  The  relative  position  of  the  bands  is  shown  schematically 
in  Fig.  4  for  a  reverse  bias  beyond  the  Kane  kink. 

IV.  INTERPRETATION 

Comparing  Figs.  1  and  2  one  sees  that  the  stress  coefficients 
r  of  Che  two  staples  differ  greatly  in  the  bias  range  of  indirect 
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tunneling,  i.e.  for  -135  <  V  <  -8  mV  and  V  >  8  mV.  We  shall 

interpret  this  difference  as  being  due  to  the  shear  induced 

shifts  of  the  (111)  valleys  and  the  nonequivalence  of  the  two  pairs 

of  valleys  for  tunneling  in  the  [110]  direction.  In  the  small  bias 

range  -8  n<V  <  V  <  +8  triV  a  very  small  current  flows.  This  must  be 

due  to  a  direct  tunneling  process  since  the  energy  difference  between 

the  Fermi  levels  on  the  n-type  and  p-type  sides  is  too  small  for 

the  emission  of  a  phonon.  Although  the  detailed  nature  of  this 

direct  tunnel  current  is  not  known,  the  relative  magnitude  of  this 

component  in  Sb,  P,  and  As-doped  germanium  tunnel  diodes  indicates^ 

that  the  origin  of  this  direct  process  is  related  to  the  impurity 

l  q 

cell  potentials  of  the  n-type  impurities.  The  fact  that  the 
stress  coefficients  of  the  two  samples  are  practically  identical 
in  this  range  indicates  furthermore  that  this  direct  tunneling 
process  cannot  be  associated  with  the  density  of  states  and  the 
band  edge  energies  of  the  individual  (111)  valleys. 

As  can  be  seen  from  Fig.  1,  both  samples  exhibit  a  sharp 
rise  in  the  stress  coefficient  at  about  140  a V  reverse  bias. 

This  effect  is  obviously  associated  with  the  onset  of  direct 
tunneling  into  the  (000)  conduction  band.  All  three  curves  are 
nearly  the  same  beyond  the  Kane  kink.  This  one  expects  since 
the  (111)  valley  contribution  to  the  current  is  negligible  in 
this  range. 
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In  the  following  we  shall  attempt  a  quantitative  comparison 
between  the  experiments  and  Kane's  theory. 

1.  Stress  Tunneling  Coefficient  Beyond  the  Kane  Kink 
One  expects  the  onset  of  direct  tunneling  to  the  (000)  conduction 
band  to  occur  at  a  reverse  bias  voltage 

Vk  "  -lEg(°0°)  "  Eg(U1)  -  £n1/e-  (5) 

Since  the  pressure  coefficient  for  the  direct  gap  E  (000)  is 

© 

0 

larger  than  that  for  the  indirect  gap  E  (111),  that  part  of  the 

© 

uniaxial  compression  X  which  corresponds  to  a  hydrostatic  pressure 
p  -  X/3  causes  the  onset  voltage  to  increase.  This  results  in 
an  anomalously  large  negative  value  of  AI/I.  As  one  goes  to  higher 
reverse  bias,  this  contribution  becomes  relatively  unimportant  and  . 
the  stress  coefficient  approaches  the  value  determined  by  the 
stress-induced  changes  of  E  (000)  and  the  combined  electron  and 

O 

hole  effective  mass  m*(000).  Figure  1  shows  that  these  changes 

are  determined  mainly  by  the  hydrostatic  pressure  part  of  the 

stress.  The  shear  part  does  not  change  E  (000).  It  deforms, 

© 

however,  the  effective  mass  spheres,  which  gives  rise  to  a  small 
contribution  which  is  different  for  the  two  diode  orientations. 

Beyond  the  Kane  kink  the  tunnel  current  is  the  sum  of  the 
direct  current  1^  and  the  Indirect  current  1^  (see  Fig.  4). 

Following  Kane  one  has 


id  -  ed  »<v  -  V  **p  (-«) 


(«) 


with 
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a  -  *d  Eg1/2(000)[«*(000j1/2/F 
and 

E  (000)  f  w  1 

D(v-vk)  -  .<v-vk)  +  -Ss-  [1-xl,tyoooy<v*vk)1j  • 

The  relative  change  of  the  tunnel  current  Is 

hi  AId  *  AI1 
"T"  Xd  +  Xl 


(7) 

(8) 

(9) 


It  can  be  seen  from  Fig.  4  that  beyond  the  Kane  kink  the  Indirect 
current  Is  an  appreciable  fraction  of  the  total  current  only  over 
a  email  bias  range.  For  the  purpose  of  explaining  the  data  we 
can  therefore  use  for  1^  and  &I^  In  Eq.  (9)  the  extrapolations  of 
the  Indirect  current  curves  into  the  bias  range  beyond  the  Kane 
kink. 


For  the  change  of  the  direct  current  one  obtains  from  Eqa.  (6) 
and  (7) »  neglecting  the  relatively  minor  change  of  C^, 


AId  "  Xd  5  AVk 


r 3  AEg(°°°)  1  Am*  (000)  AF 

-‘[2-^0001  --fJ* 


(10) 


The  first  term  in  Eq.  (10)  Is  responsible  for  the  sharp  maximum 
of  the  stress  coefficient  near  V^,  and  the  second  term  determines 
its  asymptotic  value  at  large  reverse  bias.  The  quantity  *Vkt. 


AVk  ■  *(-000  *  —  1U*  3 

where  H  000  -  12  x  lo"*2  volt  cm2/dyne  and  X  ^  ■  5  x  lo’*2 
2 

volt  cm  /dyne  are  the  pressure  coefficients  for  the  direct  and 

g 

the  indirect  band  gap,  respectively. 


(11) 
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Before  one  cen  calculate  the  theoretical  bias  dependence  of 
the  stress  coefficient  It  Is  necessary  to  estimate  a  and  the  value 
of  Am* ( 000 ) /m* ( 000 ) .  The  shear  contribution  to  the  latter  quantity 
depends  on  the  field  direction  and  hence  is  different  for  our  two 
samples.  Its  magnitude  cannot  be  obtained  without  knowing  the 
various  deformation  potentials  which  determine  the  stress-induced 
changes  of  the  effective  mass  tensors  at  the  zone  center.  Figure 
1  shows,  however,  that  the  shear  contribution  is  small.  The 
relative  change  of  the  reduced  mass  caused  by  the  hydrostatic 
pressure  can  be  estimated  from 

Am* ( 000 ) /m* { 000 )  =  AE  (000) /E  (000)  .  (12) 

8  8 

since  the  masses  involved  are  mainly  determined  by  the  interaction 

between  the  light  hold  band  and  the  conduction  band  at  the  zone 
14 

center.  The  quantity  a  was  then  determined  by  fitting  Eq.  (10) 
and  the  experimental  hydrostatic  pressure  curve  of  sanple  1  at 
large  reverse  bias.  From  this  fit,  the  value  a  ■  17.6  was 
obtained  for  V  -  -300  mV. 

As  an  independent  check,  the  bias  dependence  of  the  direct 
tunnel  current  itself  was  compared  with  Eqs.  (6)  and  (7)  using 
this  a.  The  comparison  with  experiment  is  shown  in  Fig.  5>  The 
constant  factor  of  Eq.  (6)  was  chosen  to  fit  the  magnitudes 
of  the  measured  and  the  calculated  current.  It  is  seen  that 
this  value  of  a  predicts  the  shape  of  the  I-V  characteristic 
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quite  well.  This  agreement  was  not  noted  by  Morgan  and  Kane 
because  they  assumed  the  junction  field  F  did  not  change  with 
bias.  Even  though  F  varies  by  only  about  10°/o,  the  value  of  the 
exponential  factor  varies  by  a  factor  of  6  over  the  bias  range  of 
interest. 

The  theoretical  stress  coefficient  is  compared  with  experiment 
in  Fig.  6.  It  is  seen  that  there  is  good  qualitative  agreement 
except  that  the  theoretical  maximum  is  sharper  than  the  experimental 
curve.  There  are  several  effects  which  will  cause  a  smearing  out 
of  the  theoretical  curve.  (1)  Thermal  fluctuation  will  cause 
about  a  1.5  niV  broadening.  (2)  Random  fluctuations  of  the  impurity 
concentrations  on  a  microscopic  scale  will  cause  local  fluctuations 
in  £  .  These  will  correspond  to  a  range  of  values  rather  than 
a  unique  value  as  was  assumed.  (3)  There  may  also  be  a  non-uniform 
built-in  stress  in  the  diodes.  Since  E  (000)  and  E  (111)  depend 

©  O 

differently  on  stress,  this  would  also  cause  a  spread  in  V^. 

For  pure  germanium  E  (000)  -  E  (111)  »  144  iriV.*'’  For  our 

6  © 

samples,  should  therefore  occur  at  124  vN.  The  observed 
is  clearly  larger  than  136  mV.  This  discrepancy  may  be  due  to 
a  depression  of  the  (111)  conduction  band  relative  to  the  (000) 
conduction  band  due  to  the  large  impurity  concentration^  or  it 
may  be  due  to  a  permanent  strain  at  the  junction  caused  by  the 
difference  in  lattice  constant  of  the  n-type  and  the  p-type  regions, 
or  possibly  by  the  difference  in  the  thermal  expansion  coefficient 
of  the  dot  material  and  the  germanium. 
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2.  Stress  Coefficient  in  the  Indirect  Tunneling  Range 
Except  for  the  narrow  voltage  region  at  zero  bias  (see  Fig.  2) 
the  tunneling  is  almost  entirely  indirect  for  <  V.  In  this 
bias  range  the  stress  coefficients  of  the  two  sanples  differ 
strongly.  This  difference  is  due  to  the  fact  that  in  sample  1  all 
four  valleys  are  equivalent  with  respect  to  the  junction  field 
direction,  whereas  in  sample  2  the  two  pairs  of  valleys  which  are 
shifted  with  respect  to  one  another  by  shear  have  different 
effective  mass  components  in  the  field  direction. 

The  difference  between  the  hydrostatic  and  the  uniaxial 
stress  coefficient  of  sample  1  (see  Fig.  1)  is  due  to  the  shear- 
induced  changes  of  the  effective  mass  components  in  the  field 
direction.  It  is  interesting  to  note  that  these  changes  are  of 

appreciable  magnitude.  We  find  s  5  x  10  cm  /dyne. 

A  10 

Since  m  is  the  reduced  mass  we  cannot  separate  the  contributions 
of  the  light  hole  mass  and  of  the  (111)  electron  mass.  The  fact 
that  the  shear  dependence  of  the  current  beyond  the  Kane  kink  is 
quite  small  allows  us  to  conclude  that  either  the  effective  mass 
deformation  occurs  mostly  on  the  (111)  ellipsoids  or  that  there 
is  an  accidental  cancellation  of  the  deformations  of  the  light 
holes  and  the  (000)  conduction  band  extremum. 

The  bias  dependence  of  the  hydrostatic  pressure  coefficient 
shows  a  clearly  resolved  structure  at  small  bias  voltages  which 
is  symmetric  with  respect  to  the  forward  and  reverse  bias  direction. 


-  19  - 


The  value  of  the  pressure  coefficient  ■  AI/3lp  is  about  « 

-11  2, 

-7.5  x  10  cm  /dyne  in  the  small  bias  range  of  direct  tunneling. 

Tp  changes  to  -4  x  10  cm  /dyne  at  V  =+  6  mV  before  it  reaches 
a  constant  value  of  -5.5  x  10  Near  V  =  +  30  rxN ,  which  is  the 

onset  of  the  LA  phonon,  v  changes  again  to  about  *  -6.5  x  lo"^ 
cm  /dyne.  This  structure  is  also  indicated  although  not  as  clearly 
by  the  uniaxial  stress  coefficient  of  sample  1.  It  is  to  be  expected 
that  the  stress  coefficient  of  the  direct  tunneling  current  at 
very  small  biases  is  different  from  that  of  the  Indirect  tunneling 
current.  We  cannot  explain,  however,  the  remaining  part  of  the 
structure. 

The  pressure  coefficient  of  sample  1  in  the  Indirect  tunneling 
range  is  due  to  the  change  of  the  indirect  band  gap  E  (111)  and 

O 

4* 

that  of  the  reduced  mass  m  .  Neglecting  the  minor  stress  variations 
of  and  D  one  can  write  from  Eq.  (2) 


3  AV1U)  ,  1  Am4  AF 
l2  E*(lli)  +  2  '  F 


m 


(13) 


One  may  estimate  0  from  Eq.  (13)  and  the  measured  w  « 

6  x  10  cm  /dyne  of  sample  1  in  the  forward  bias  range.  Here 
again  it  is  difficult  to  estimate  Am4/m4.  The  reduced  mass  m4 
in  Eq.  (13)  contains  the  light  hole  mass  and  the  component  of  the 
(111)  valley  mass  tensor  in  the  field  direction.  The  change  of 
the  light  hole  mess  can  be  estimated  from  Eq.  (12)  where 
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.  dE  (000)  ,2  2 

^  ^QQ0'j  — "jp «  12.4  x  10  cm  /dyne.  The  change  of  Che 

transverse  electron  mass,  which  predominantly  determines  the 

component  of  the  electron  mass  in  the  field  direction,  is  related 


.-12 


to  the  band  gap  change  at  the  zone  face  at  [111]  which  is 
dE  A 

2 - -  4  x  10~*2  cm2/dyne.^  Hence  *  (8  +  5)  x  10 

go  ® 

cm  /dyne  may  be  a  reasonable  estimate.  Using  this  value  and  the 
measured  pressure  coefficient  one  obtains  0  «  16  +  3  for  V  -  +60  vtt . 
The  same  calculation  yields  the  value  0-20  +  4  for  V  «.  -70  vN. 

This  increase  of  0  with  decreasing  bias  voltage  disagrees  strongly 
with  Eqs.  (2)  and  (3)  which  predict  instead  a  slight  decrease  of 
S.18 


One  may  think  that  this  discrepancy  is  due  to  some  simplifying 

assusptlons  on  which  the  derivation  of  Eqs.  (2)  and  (3)  is  based, 

in  particular  the  assumptions  of  (1)  a  constant  junction  field, 

(ii)  a  parabolic  shape  of  the  energy  bands  in  k-space,  and  (ill) 

19 

an  abrupt  transition  from  n-type  to  p-type  doping.  Nathan, 
however,  has  calculated  the  tunneling  exponent  0  using  the  more 
realistic  field  of  a  graded  junction  and  including  non-parabolic 
effects.  He  finds  that  both  the  variation  of  the  field  with 
position  in  the  Junction  and  a  graded  impurity  distribution,  causa 
0  to  decrease  with  decreasing  bias  voltage  even  faster  than 
predicted  by  Eqs.  (2)  and  (3)  and  that  the  nonparabolic  effects 
do  not  Influence  the  bias  dependence  of  0  appreciably.  In  addition 


l 
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to  this  Nathan  finds  that  for  interpreting  the  current -voltage 
characteristic  of  germanium  diodes  at  297°K  one  requires  p  to 
decrease  with  decreasing  voltage.  In  comparing  the  current-voltage 
characteristic  of  the  direct  tunneling  process  beyond  the  Kane 
kink  with  Eqs.  (2)  and  (3)  in  Fig.  5  we  also  concluded  that  a 
decreases  with  decreasing  bias. 

In  view  of  this  evidence  it  is  difficult  to  understand  the 
increase  of  0  with  decreasing  bias  voltage  as  observed  in  the 
hydrostatic  pressure  experiments. 

Since  one  is  dealing  with  the  same  diode  one  should  be  able 
to  calculate  0  from  the  previously  determined  a  =  17.6  at  V  =  -300  ntf. 
Using  Eqs.  (1),  (2),  and  Kane's  value  X^Ad  =  one  obtains  for 
our  case  p/a  «  1.5  and  hence  p  *  26.5  at  V  *  -300  mV.  This  value 
of  p  compares  well  with  P  «  23.5  +  5  obtained  by  extrapolating 
the  observed  bias  dependence  of  p  to  V  *  -300  vN.  However,  in 
view  of  the  unexpected  bias  dependence  of  P,  this  agreement  may 
be  fortuitous. 

The  contribution  to  the  stress  coefficient  of  satqple  2  which 
arises  from  the  non-equivalence  of  the  two  pairs  of  valleys 
consists  of  two  parts.  (1)  The  shear  causes  some  electrons  to 
transfer  from  the  two  valleys  which  are  raised  in  energy  into  the 
two  which  are  lowered.  Those  which  are  lowered  are  the  easy 
valleys,  i.e.  they  have  the  smaller  a*  in  the  tunneling  direction. 


r 
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Hence  we  expect  thl*  pert  to  give  e  positive  contribution  to  ir. 

(2)  The  shift  of  the  valleys  will  change  also  the  tunneling  prob¬ 
ability  because  the  energy  gap  for  a  particular  valley  is  changed. 
Since  the  easy  valleys  are  lowered  in  energy,  their  energy  gap 
decreases  under  shear.  This  again  results  in  a  positive  contribution 
to  r. 


Since  the  first  effect  involves  the  density  of  states  factor 
D  one  has  to  write  the  current  as  a  sum  over  the  j  valley 
contributions  and  the  k  phonon  contributions 

xJk 

where 


xjk  *  \  °(ev  Ve)  exp  (‘V  (l5) 

with 

Bj  "  \  Eg3/2(U1)  ®j*1/2/F  (16) 


The  plus  and  minus  signs  in  Eq.  (15)  indicate  that  the  phonon 
voltage  ^ko^/e  has  to  be  subtracted  from  the  magnitudes  of  both 
the  forward  and  the  reverse  bias  voltage  for  calculating  D^.  The 

coefficients  take  into  account  the  different  electron-phonon 

12 

Interaction  strengths  for  the  two  phonons.  From  Fig.  3  we 
determined  Ag/A^  -  2.23  where  Ag  refers  to  the  higher  energy 
phonon. 

The  stress  coefficient  due  to  the  non-equivalence  of  the 
velleys  will  then  be 

anx  ■  2^  *V2jk  V  x  <17> 


V 
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From  Eq.  (13)  one  obtains: 

Mjk  -  \  «p  (-Bj)  (^n)j  -  |  AEgj(UlJ  (18) 

Since  we  consider  In  (17)  the  pure  shear  part  of  the  stress 
only  we  have 

A  Eg  (111)  -  0  (19) 

j  J 

and  tone* 

(Af_)  -  -A  E  (111)  .  (20) 

J  si 

For  uniaxial  compression  along  [lTo]  one  finds 

A  E  (111)  -  -A  Eg^(lll)  -  ^  E2  Sw  X  (21) 

Here  the  subscript  1  refers  to  the  easy  valleys  j  =  1,  2  and 

subscript  3  to  the  hard  valleys  j  ■  3,  4.  E2  =  19  eV  is  the 

20  -12 

deformation  potential  for  pure  shear  and  =  1.47  x  10 

2  21 
cm  /dyne  Is  the  elastic  compliance  constant. 

In  the  second  term  of  Eq.  (18),  which  is  due  to  the  change 

in  tunneling  probability,  we  have  neglected  the  shear- induced  change 

of  the  valley  masses.  Since  the  reduced  mass  is  determined  more 

by  the  light  hole  mass  than  by  the  electron  mass,  this  effect  is 

not  large.  Any  change  in  the  light  hole  mass,  however,  will 

affect  both  samples  in  the  same  way  and  hence  will  not  contribute 

to  the  difference  of  r  of  the  two  senples. 


w 
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Because  of  the  uncertainties  involved  in  the  previous 
determinations  of  a  and  {3,  we  have  adjusted  the  value  of  (3^  so 
that  the  magnitude  of  Eq.  (17)  agrees  with  the  difference  in 

of  the  two  samples  in  the  region  where  only  the  first  phonon 

a  a  1/2 

contributes.  With  =  (m  ^/m  )  '  we  obtained  p^  ■  16  and 

■  18,  both  values  quoted  for  V  =  0. 

Figure  7  shows  the  comparison  of  Eq.  (17)  with  the  difference 
in  ir  of  the  two  samples.  Although  the  structure  of  the  bias 
dependence  is  reproduced  quite  well  by  the  theoretical  curve, 
the  agreement  is  only  qualitative.  The  effect  of  the  shift  in 
valley  population  can  clearly  be  seen  in  the  region  between  .008 
and  .030  volts  bias.  The  discrepancy  becomes  large,  however,  as 
soon  as  the  higher  energy  phonon  starts  to  contribute  to  the 
tunneling  current.  Here  the  experimental  data  falls  off  abruptly 
from  the  theoretical  curve.  There  is  no  known  reason  why  the 
two  phonons  should  behave  differently,  but  in  order  to  see  what 
would  happen,  a  curve  was  plotted  as  though  the  valleys  were 
equivalent  for  tunneling  involving  the  higher  energy  phonons. 

This  curve  is  shown  dashed  in  Fig.  7-  As  can  be  seen  from  the 
forward  bias  data,  this  modification  is  too  drastic.  The  data 
in  the  reverse  bias  direction  is  in  better  agreement,  but  it  is 
hard  to  imagine  any  mechanism  which  would  distinguish  forward 
from  reverse  bias.  The  fall-off  of  the  data  in  the  reverse  bias 
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direction  will  be  affected  by  exactly  the  same  smearing  out 
mechanisms  which  are  active  beyond  the  Kane  kink.  If  these 
mechanisms  are  indeed  responsible  for  the  slow  fall-off  in  that 
region,  similar  behavior  in  this  bias  region  is  to  be  expected 
also.  This  smearing  would  not,  however,  account  for  the  fact  that 
the  observed  difference  between  the  ir  of  the  two  samples  near  V  - 
-0.13  volts  is  appreciably  smaller  than  the  calculated  shear  co¬ 
efficient.  We  therefore  conclude  that  the  variation  of  the 
inherent  tunneling  probability  per  electron  is  not  as  sensitive  to 
shear  as  would  be  expected.  Within  the  framework  of  the  theory 
the  only  possibilities  which  can  account  for  this  are:  (a)  that 
the  chosen  value  for  0  is  too  low;  (b)  that  the  effect  of  shear 
on  the  reduced  effective  mass  is  quite  large  and  in  the  opposite 
direction  as  the  effect  of  the  valley  energy  shifts. 

There  is  no  value  of  p  which  is  simultaneously  consistent 
with  the  shear  data  and  the  hydrostatic  pressure  data.  In  view 
of  the  general  agreement  between  the  various  Independent  calculations 
of  P,  possibility  (a)  is  unlikely.  Possibility  (b)  can  be  eliminated 
also  because  of  the  similarity  of  ir  for  sample  1  and  t^.  Further¬ 
more,  the  mechanism  responsible  for  the  discrepancy  beyond  30  n tJ 
in  the  forward  bias  region  is  undoubtedly  also  operating  in  this 
bias  region.  It  therefore  seems  that  a  different  mechanism  not 
accounted  for  by  the  present  theory  contributes  to  the  shear 
•tress  coefficient. 
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3.  Comparison  of  Absolut*  Magnitudes  of  Tunneling  Parameters 
with  Theory 

The  physical  properties  of  a  diode  determine  the  stress 
coefficients  through  the  parameters  a  and  0.  These  parameters 
could  not  be  determined  unambiguously  from  the  stress  coefficients 
because  the  stress-induced  changes  of  the  relevant  effective  masses 
are  not  known.  These  measurements  would  yield  these  effective 
mass  changes  if  it  were  possible  to  determine  the  parameters  a  and 
0  by  some  other  method.  The  theory  leading  to  Eqs.  (1)  and  (2) 
does  not  determine  these  parameters  accurately  enough  since  it 
is  based  on  the  assumption  of  a  constant  junction  field  and 
neglects  the  spatial  distribution  and  fluctuations  of  the  impurity 
concentrations  near  the  junction. 

It  is  nevertheless  interesting  to  compare  and  evaluate  the 
values  of  a  and  0  predicted  by  the  constant  junction  field  theory 
with  the  values  obtained  above  from  the  stress  coefficients  using 
reasonable  estimates  for  Am  /m  .  The  coefficients  a  and  0  can 
be  obtained  from  the  theory,  (1)  by  direct  computation  using 
Eqs.  (1),  (2),  and  (3),  and  (ii)  by  calculating  C  and  D  and  using 
the  measured  I  for  a  given  bias  voltage.  This  was  done  with  the 
following  choice  of  parameters  appropriate  for  our  diodes, 

£p  -  0.15  eV,  -  0.020  eV,  p  -  6  x  1019,  n-5-5*  1018, 
junction  area  -  0.002  cm2,  *  -  16,  Id  ■  0.370  A  at  V  >  -300  stf. 
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The  values  for  a  and  $  obtained  by  method  (1)  are  Hated  in 
the  first  row  of  Table  I.  The  third  row  gives  a  as  calculated  by 
method  (ii).  The  value  of  0  could  not  be  obtained  in  this  way 
because  the  electron  phonon  coupling  constants  are  not  known.  The 
fourth  and  fifth  rows  list  the  values  obtained  from  the  experimentally 
observed  hydrostatic  pressure  coefficients  as  explained  in  the 
previous  section.  The  last  row  gives  the  value  of  0  used  to  fit 
the  shear  contribution  to  the  stress  coefficient  of  sample  2  in 
the  bias  range  of  indirect  tunneling.  The  bias  voltages  at  which 
the  listed  quantities  were  obtained  are  included  in  Table  I. 

It  is  seen  that  a  and  3  as  computed  from  the  constant  field 
expressions  Eq.  (1),  (2),  and  (3)  are  appreciably  smaller  than 
the  other  values.  This  is  to  be  expected  since  a  real  diode  will 
have  a  more  diffused  distribution  of  impurities  and  hence  a  smaller 

impurity  concentration  near  the  junction  than  in  the  bulk  material. 

22 

This  yields  a  wider  junction.  Following  Meyerhofer  et.  al.  we 
assume  the  donor  concentration  near  the  junction  to  be  half  of  that 

of  the  bulk  and  obtain  by  direct  computation  of  a  and  0  the  values 

23 

quoted  in  the  second  row  of  Table  I. 

The  uncertainty  of  the  experimental  values  stems  from  the 
difficulty  of  estimating  Am* /m*  and  also  from  the  fact  that  for 
the  pressure  coefficients  of  the  energy  gaps  the  values  measured 
on  pure  Mterial  at  300°K  were  used.  It  is  hoped  to  obtain  the 
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temperature  dependence  of  these  pressure  coefficients  by  studying 
the  stress  tunneling  coefficients  at  elevated  temperatures. 

V.  SUMMARY  AND  CONCLUSION 

The  magnitude  and  the  orientation  dependence  of  the  effect  of 
stress  on  the  tunneling  current  enables  one  to  distinguish  three 
different  tunneling  processes  in  Sb-doped  germanium. 

(1)  In  the  small  bias  range  -8  niV  <  V  <  +8  mV  the  tunneling 
is  direct.  The  detailed  nature  of  the  tunneling  process  in  this 
range  is  not  known.  The  current  is  not  affected  by  the  relative 
shifts  of  the  (111)  valleys  under  shear  and  hence  cannot  arise 
from  a  sum  of  independent  (111)  valley  contributions. 

(2)  Beyond  about  8  nff  at  helium  temperatures,  but  before 
the  onset  of  direct  tunneling  into  the  (000)  conduction  band,  the 
tunneling  process  is  phonon  assisted.  Here  a  large  additional 
contribution  to  the  stress  coefficient  is  observed  for  a  shear 
stress  which  lifts  the  degeneracy  of  the  (111)  valleys  when  the 
orientation  of  the  junction  field  is  such  that  the  effective 

mass  components  of  the  valleys  in  the  field  direction  are  different. 
This  demonstrates  clearly  the  anisotropy  of  the  tunneling  probability 
when  the  effective  mass  is  anisotropic  as  predicted  by  the  theory. 

(3)  For  biases  V  <  -140  vN  direct  tunneling  into  the  (000) 
conduction  band  is  observed.  Near  the  onset  voltage  the 
magnitude  of  the  stress  coefficients  lncreeaes  sharply  because 
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of  the  pressure  dependence  of  V^*  This  sharp  increase  allows  an 

accurate  determination  of  E  (000)  -  E  (111),  the  energy  separation 

8  8 

between  the  conduction  band  edges  at  (000)  and  (111),  respectively. 
We  find  for  this  the  value  0.160  +  0.005  eV  which  is  about  0.015  eV 
larger  than  that  for  pure  material.  This  larger  value  of  the 
conduction  band  separation  agrees  with  recent  results  of  infrared 
absorption  measurements^  in  degenerate  germanium.  The  sharp  rise 
of  the  stress  coefficients  near  V  =  indicates  that  the  (000) 
conduction  band  edge  is  quite  well  defined  despite  the  large 
impurity  concentration. 

The  major  features  of  the  bias  and  orientation  dependencies 
of  the  stress  coefficients  agree  with  the  present  theory  of 
tunneling.  There  are,  however,  several  interesting  observations 
which  remain  unexplained  at  present.  One  is  the  structure  in  the 
bias  dependence  of  the  pressure  coefficient  irp.  The  magnitude 
of  Tp  changes  near  V  =  +  8  vN  and  V  =  +  30  vtN.  These  voltages 
agree  closely  with  the  onsets  of  the  TA  and  the  LA  phonon 
contributions.  The  second  is  the  disagreement  between  the 
theoretically  predicted  shear  contribution  to  ir  with  the  experiments 
in  the  bias  range  where  LA  phonons  can  be  emitted.  The  third 
observation  which  needs  some  further  study  is  the  bias  dependence 
of  the  indirect  tunneling  exponent  $  which  is  found  to  be  opposite 
to  the  prediction  of  theory. 
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TABLE  I 

Values  of  Tunneling  Exponents  a  and  0 


a 

Bias  (ntf) 

P 

Bias  (nf 

Theory 

8.5 

-300 

12.5 

-300 

Theory  (see  text) 

12 

-300 

17-5 

-300 

Theory  and  1^ 

12.5+1 

-300 

Pressure  Exp. 

17.6 

-300 

16  +  3 

+  6o 

Pressure  Exp. 

20+4 

-  70 

Shear  Exp. 

17  +  1 

0 
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FIGURE  CAPTIONS 


Fit*  1 


Fit*  2 


FI*.  3 


FI*.  4 


FI*. 


FI*.  6 


FI*.  7 


Stress  coefficients  for  uniexiel  Elio]  compression  end 
hydrostatic  pressure  as  a  function  of  bias  volta*a  at 
4.2°K. 

Current -vo It age  characteristic  of  one  diode  of  sample  2 
at  4.2°K.  Note  the  onsets  of  the  TA  and  LA  phonon  con¬ 
tributions  near  the  bias  voltages  +  8  VN  and  +  28  wN, 
respectively. 

Current  as  a  function  of  reverse  bias  of  one  diode  of 
sample  2  at  4.2°K.  Note  the  sharp  increase  in  current 
near  V  -  -140  mV. 

Relative  position  of  the  bands  and  Fermi  levels  in  a 
germanium  tunnel  diode  for  a  reverse  bias  beyond  the 
Kane  kink. 

Comparison  between  theory  and  experiment  of  the  bias 
dependence  of  the  direct  tunnel  current  beyond  the  Kane 
kink.  The  factor  C^  of  Eq.  (1)  (see  text)  has  been 
chosen  to  fit  the  absolute  magnitude  of  the  theoretical 
and  experisiental  1^  near  -220  vN  reverse  bias. 
Comparison  between  theoretical  and  experimental  stress 
coefficient  for  reverse  bias  voltages  beyond  the  Kane 
kink. 

The  difference  between  the  uniaxial  stress  coefficients 
of  samples  1  and  2.  Comparison  between  theory  and 
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•xperiment.  The  heavy  curve  represents  the  theory 
based  on  the  anisotropy  of  the  tunneling  probability 
in  the  whole  bias  range  of  indirect  tunneling.  The 
dotted  curve  was  calculated  assuming  that  the  (111) 
valleys  are  equivalent  for  tunneling  involving  the 
higher  energy  phonon. 
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B.  Carrier  Mobility  and  Shallow  Impurity  States  in  ZnSe  and  Zaire 
(M.  Atrea  and  B.  Segall) 

I .  INTRODUCTION 

Several  of  the  II-VI  compounds  have  been  under  investigation 
by  luminescence  and  photoconduct ive  techniques  for  a  great  many 
years.  Much  of  this  work  was  directed  toward  the  goal  of  under¬ 
standing  the  role  of  the  impurities  and  of  identifying  the  states 
associated  with  them.  Despite  this  work,  however,  the  identification 
of  the  specific  Impurity  states  is  still  a  matter  of  considerable 
debate  and  uncertainty.  At  present  there  is  still  little  conclusive 
evidence  about  the  identity  of  any  shallow  impuritystates,  which, 
for  example,  play  an  important  role  in  the  optical  properties  of 
these  crystals. 

The  primary  aim  of  the  present  work  was  to  obtain  information 
about  the  shallow  donor  and  acceptor  levels  in  some  of  the  II-VI 
compounds  by  means  of  electrical  transport  studies.  In  the  past 
extensive  electrical  measurements  have  not  generally  been  possible 
due  to  the  unavailability  of  suitable  material.  However,  considerable 
advances  have  recently  been  made  in  the  preparation  of  some  of 

these  materials,  including  the  compounds  that  are  the  subject  of 

12; 

the  present  study,  ZnSe  and  ZnTe.  '  It  is,  of  course,  well 
known  that  in  addition  to  yielding  information  about  the  energy 
levels  of  the  defects,  the  electrical  transport  measurements  provide 
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a  valuable  means  for  studying  the  mechanism  by  which  the  carriers 
are  scattered.  This  question  was  considered  In  some  detail  for 
the  two  compounds  studied. 

In  earlier  work,  the  results  of  optical  and  preliminary 

4  5 

electrical  measurements  on  ZnSe  were  reported.  The  electrical 
measurements  of  the  n-type  material  have  been  extended  in  the 

present  work.  Besides  the  previously  described  material  with 

5  6 

0.19  ev  deep  donor  levels,  and  the  more  strongly  n-type  material 

15 

containing  donors  with  ionization  energies  of  0.01  eV,  it  has 
now  been  possible  to  prepare  ZnSe  exhibiting  degenerate  conduction 
characteristics . 

We  have  also  extended  the  investigation  of  the  transport 
properties  of  p-type  ZnTe.  The  study  of  this  compound  is  of 
considerable  interest  since  the  knowledge  gained  about  the  transport 
of  holes  in  the  valence  bands  and  of  the  acceptor  states  would  be 
expected  to  be  helpful  in  understanding  the  corresponding  properties 
and  states  of  the  other  II -VI  semiconductors.  In  regard  to  the  study 
of  acceptor  levels,  significant  progress  was  made  possible  as  a 
result  of  the  successful  purification  by  the  zinc  extraction 
technique  recently  described  by  Aven  and  Woodbury.  By  impurity 
diffusion  and  extraction  cycles,  in  conjunction  with  radioactive 
tracer  studies,  it  was  possible  to  identify  the  acceptor  levels 
associated  with  the  most  prominent  acceptor  lmpurlties--Cu,  Ag, 
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and  Au.  Concomitantly,  the  acceptor  state  of  the  native  defect, 
which  we  believe  to  be  the  zinc  vacancy,  could  be  Identified.  It 
Is  to  be  emphasized  that  analogous  states  are  believed  to  be 
associated  with  the  other  crystals  of  this  family  of  semiconductors, 
but  until  this  time  the  identification  of  these  states  has  not  been 
conclusively  established. 

From  the  study  of  the  Hall  mobility,  the  intrinsic  (or  lattice) 
mobility  for  the  two  materials  could  be  determined  for  temperatures 
above  about  100°K.  The  various  possible  lattice  scattering 
mechanisms  were  considered  with  the  aim  of  determining  which  of 
them  play  the  major  role  in  limiting  the  intrinsic  mobility.  For 
n-type  ZnSe  it  is  found  that  the  electrons  are  scattered  such  more 
strongly  by  the  polar  scattering  of  optical  phonons  than  by  other 
mechanisms.  A  similar  conclusion  is  probably  also  true  for  ZnTe, 
unless  the  parameter  characterizing  the  strength  of  the  non -polar 
coupling  to  the  optical  mode  phonons  is  much  larger  than  presently 
believed. 

II.  EXPERIMENTAL 

The  ZnSe  and  ZnTe  crystals  were  grown  by  the  vapor  growth 

g 

technique  first  described  by  Greene  et  al  and  modified  by  PJtper 

Q 

and  Polich.  The  starting  material  for  ZnSe  was  General  Electric 
Chemical  Products  Plant  ZnSe  powder  which  was  first  sintered  in 
the  presence  of  excess  Se  in  hydrogen  atmosphere.  ZnTe  waa 


r 
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synthesized  by  heating  together  high  purity  Zn  obtained  from  the 
Consolidated  Mining  and  Smelting  Co.  and  high  purity  Te  obtained 
from  the  American  Smelting  and  Refining  Co.  The  single  crystals 
were  grown  by  passing  the  sintered  charge  at  the  rate  of  approxi¬ 
mately  0.2  iwn/hr  through  a  furnace  heated  to  1350°C  for  ZnSe  and 

to  1050°C  for  ZnTe.  The  product  was  usually  a  boule  of  several 

•a  3 

cm3  containing  single  crystal  regions  of  a  tenth  to  one  cm  . 

The  majority  of  the  crystals  were  subjected  to  the  purifi¬ 
cation  technique  in  which  the  crystals  are  heated  in  contact  with 
molten  Zn  at  900-l000°C  for  one  or  two  days.  This  treatment  has 
been  shown^  (by  using  radioactive  tracer  techniques)  to  reduce 
the  concentration  of  the  acceptor  impurities  Cu  and  Ag  in  II-VI 

14  ic  _o 

compounds  to  the  10  to  10  J  cm  range.  With  the  exception  of 
Cl-doped  ZnSe,  all  the  doped  crystals  were  oH-ained  by  diffusing 
impurities  into  crystal  purified  by  the  sol^nt  extraction  method. 
Chlorine-doped  ZnSe  was  prepared  by  adding  SrClg  to  the  starting 
material  prior  to  the  crystal  growth. 

The  Hall  measurements  were  performed  with  sample  bars  8-10  mm 

p 

long  and  1-4  mn  in  cross  sectional  area,  cut  from  the  boules 
described  above.  Prior  to  electroding  the  bars  were  chemically 
polished  in  hot  concentrated  NaOH.  The  electrodes  to  ZnSe  were 
In  dots  fused  in  under  hydrogen  atmosphere  at  300°C.  The  electrodes 
to  ZnTe  were  chemically  formed  Au  patches,  to  which  Ft  wires  were 
soldered  with  In-Ag  alloy. 
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The  Hell  effect  apparatus  used  for  electrical  measurements 
covered  the  range  from  liquid  hydrogen  temperature  to  approxi¬ 
mately  150°C.  The  measurements  were  performed  by  the  standard 
six  electrode  technique  widely  used  on  semiconductor  materials. 

The  voltage  drop  across  the  crystal  was  usually  of  the  order  of 
0.1  V.  The  magnetic  fields  used  were  between  7  and  8  kgauss. 

To  assure  that  no  irreversible  changes  had  occurred  in  the  bulk 
of  the  crystals  or  the  electrodes  during  the  thermal  cycling, 
several  of  the  crystals  were  remeasured  over  the  entire  covered 
temperature  range.  A  few  crystals  were  removed  from  the  apparatus, 
re-electroded  and  remeasured.  In  no  cases  were  any  irreversible 
effects  observed.  However,  with  all  crystals  a  minimum  temperature 
was  always  encountered  below  which  the  measurements  became  unreliable 
because  of  the  onaet  of  non-ohmic  behavior  of  the  electrical 
contacts.  Thus  for  any  given  crystal  the  temperature  range  for 
which  the  measurements  are  given  was  determined  by  the  lowest 
temperature  at  which  the  Hall  voltage  and  the  potential  drop  across 
the  crystal  were  linear  functions  of  the  current  through  the  crystal. 

111.  RESULTS  AND  DISCUSSION 

Before  discussing  and  analyzing  the  experimental  results,  it 
is  useful  to  state  the  assumptions  about  the  band  structure  of 
ZnSe  and  ZnTe  that  will  be  used.  It  appears  that  both  of  these 
compounds  are  direct  transition  semiconductors,  presumably  having 
the  absolute  extreme  of  their  conduction  end  valence  bands  at  k  •  0. 
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These  conclusions  are  based  primarily  on  the  absorption  spectra  * 

(particularly  on  the  sharp  exciton  lines  which  are  not  appreciably 

11  12 

broadened  by  autoionization)  and  on  the  emission  spectra  ’  which 

13 

exhibit  sharp  lines  very  close  to  the  band  edges.  Since  we  will 
only  be  concerned  with  n-type  ZnSe  and  p-type  ZnTe,  only  the 
assumptions  about  the  location  of  the  conduction  band  minimum  of 
the  former  and  the  valence  band  maximum  of  the  latter  will  have  a 
bearing  on  the  following  discussion.  The  chief  effect  of  this 
assumption  will  be  to  exclude  the  possibility  of  intervalley 
scattering. 

In  addition,  it  will  be  necessary  to  assume  values  for  the 
effective  masses  for  the  carriers.  On  the  basis  of  a  reduced 

exciton  mass  of  (0.10  +  .03)  m  as  deduced  in  the  optical  studies 

4  * 

of  ZnSe  ,  we  take  m  -  0.15  m  as  a  reasonable  estimate  of  the 

electron  mass.  Since  Information  relating  to  the  hole  mass  in 

ZnTe  is  unavailable,  we  are  forced  to  make  a  more  arbitrary 

14 

choice  in  this  case  based  on  the  known  hole  masses  in  CdS  and 
15 

CdSe.  We  take  the  value  of  0.6  m  for  the  density  of  states 
mass.  Inasmuch  as  we  are  not  seeking  to  precisely  fit  the  data, 
this  estimate  of  the  hole  mass  should  be  adequate  for  our  purpose. 

In  analyzing  the  carrier  concentration  data  from  a  single 
level,  we  will  use  the  wall-known  formula1**  for  non-degenerate 
statistics,  e.g. 

P(P  +  V  Nv 

dt.-V*5  " "»  ,xp(’E«/kI)  (1) 
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£or  the  hole  concentration  p.  Here  Nv  =  2(2rmvkTh  )  . 

being  the  density  of  states  mass  and  g  the  degeneracy  factor 

which  depends  on  the  nature  of  the  impurity  state  and  the  band 

edge  involved.  Later  when  considering  the  carrier  concentration 

associated  with  an  imperfection  which  we  believe  to  be  a  double 

acceptor,  we  will  set  up  the  generalization  of  (1)  appropriate 

to  that  case.  In  the  following,  we  will  take  the  carrier  con* 

cent  rat  ions  p  and  n  to  be  given  by  j^e)  1  for  simplicity,  where 

Rjl  is  the  Hall  constant  and  e  the  electronic  charge.  The  factor 

r  ■  u„/u.j,^  the  ratio  of  Hall  to  drift  mobilities  is  thus  taken 
H  a 

to  be  unity.  This  is,  in  fact,  a  fairly  good  approximation  since 

17 

r  is  reasonably  close  to  unity  for  polar  scattering  which  we 
will  show  below  to  be  the  principal  scattering  mechanisai  in  the 
temperature  range  of  primary  interest  here. 

1.  ZnTe  -  Impurity  and  Native  Defect  Levels 
In  the  following  we  will  discuss  the  results  of  Hall 
measurements  on  several  ZnTe  crystals  prepared  in  such  a  way  as 
to  produce  a  dominance  of  one  type  of  impurity  (including  native 
defect)  over  all  others.  With  such  crystals  it  is  a  straight* 
forward  matter  to  identify  the  various  defects  giving  rise  to  the 
energy  levels  found  by  the  electrical  measurements. 

A  survey  of  the  behavior  of  the  Hall  constant  vs.  T*  for 
typical  undoped  crystals  and  crystals  doped  with  the  noble  metals 
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Is  provided  by  Fig.  1.  From  the  curve  for  the  undoped  sample ,  the 
p-type  nature  of  which  we  believe  to  be  due  to  a  Zn  vacancy,  an 
Ionization  energy  of  0.048  eV  is  found.  The  acceptor  ionization 
energies  for  the  Cu,  Ag,  and  Au-doped  samples  are  found  to  be 
0.15  eV,  0.11  eV,  and  0.22  eV  respectively. 

Figure  2  shows  more  detailed  data  for  a  crystal  which  was 
found  by  chemical  analysis  to  contain  Cu  at  a  concentration  of 

16  -  o 

(2  +  1)  x  10  cm  .  The  behavior  at  low  T  corresponds  to  the 
freeze-out  of  the  holes  into  an  acceptor  level  at  0.15  eV  as  in 
the  Cu  curve  of  Fig.  1.  Using  the  value  of  g  *  4  which  is  appropriate 
for  a  simple  acceptor  associated  with  the  four-fold  degenerate 

ifl 

valence  band  maximum  it  is  possible  to  satisfactorily  fit  the 
data  with  »  3.4  x  10*^  cm’^  and  «  3.1  x  lO1^  cm  In 
Fig.  2  the  experimental  points  are  indicated  by  circles,  and  the 
theoretical  values  by  the  solid  curve.  It  should  be  noted  that 
the  above  value  of  NA  agrees  reasonably  well  with  the  result  of 
the  chemical  analysis. 

Figure  3  shows  the  free  hole  concentration  vs  T**  for  four 
undoped  ZnTe  crystals  fired  under  different  Te  pressures.  The 
shapes  of  these  curves  are  sufficiently  similar  to  indicate  that 
the  same  level  is  involved  in  all  the  crystals.  The  energy  value 
that  corresponds  to  the  low  temperature  end  of  the  curves  is 
0.048  eV  (the  temperature  range  for  the  two  lower  curves  being 
admittedly  very  small) .  The  vertical  displacement  of  theee  curvee 


t 
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from  one  another  indicates  that  the  free  hole  concentration  at  any 
temperature  increases  with  increasing  partial  pressure  of  Te.  The 
Increasing  partial  pressure  of  Te  increases  the  concentration  of 
Zn  vacancies,  which  in  turn  increases  the  free  hole  concentration* 
These  results  are  consistent  with  our  belief  that  the  native 
acceptor  center  is  the  Zn  vacancy. 

The  above  evidence  per  se  does  not  preclude  the  possibility 
that  the  acceptor  is  interstitial  Te.  However,  this  possibility 
can  probably  be  ruled  out  on  the  basis  of  the  very  large  size  of 
the  Te  atom  and  also  because  it  appears  unlikely  that  interstitial 
Te  could  bind  extra  electrons.  The  tendency  of  ZnTe  to  form  an 
excess  of  Zn  over  Te  vacancies,  on  the  other  hand,  is  consistent 
with  the  greater  volatility  of  Zn. 

The  Zn  vacancy  is  expected  to  be  a  double  acceptor.  In 
contrast  to  the  curves  for  the  Cu,  Ag,  and  Au-doped  samples,  the 
curves  for  the  undoped  samples  have  a  shape  above  about  200°K 
which  suggests  the  ionization  of  holes  from  a  deeper  level.  In 
order  to  check  the  identification  and  determine  the  ionization 
energy  of  the  deeper  level,  it  would  be  desirable  to  have  hole 
concentration  data  available  over  a  wide  high  temperature  range, 
and  to  fit  these  data  with  the  expression  for  p  appropriate  to 
the  case  of  a  double  acceptor  associated  with  the  fourfold 
degenerate  valence  band.  Although  the  present  data  is  perhaps 
not  as  extensive  as  would  bs  desired,  it  still  is  useful  for  a 
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preliminary  investigation  of  this  point.  To  obtain  the  equation 
for  p  we  use  for  the  probabilities  that  the  double  acceptor  ia  in 

the  neutral,  single,  or  doubly  ionised  state  the  expressions  derived 

18  I 

by  Teltler  and  Wallis,  which  are 

Na°/Na  -  6  R0_1  ^exP[2a  +  &(Ei  +  ] 

Na~/Na  "  4  R0*1  +  p  E2J  (2) 

Nfl2_/N  -  R- 1 
a  a  o 

with  N  °  +  N  ~  +  N  2"  -  N  ,  a  -  -E_/kT,  and  £  *  1/kT.  For  non- 
a  a  a  a  f 

degenerate  statistics,  the  charge  neutrality  condition  p  +  Nd  » 

.  2- 

+  2Na  leads  to  the  following  cubic  equation  for  p: 

6p3  +  p2  [4n^1}  +  6Nd]  +  p  [N^1}n|2)  +  4N^L)  (N^)] 

+  l)  n|2)  (Nd-2Nm)  -  0  (3) 

where  -  Ny  exp  (-E^/KT)  for  i  ■  1  and  2.  At  temperatures  low 

(2) 

enough  so  that  is  negligible  Eq.  (3)  reduces  to  Eq.  (1)  with 

iq 

g  »  3/2  which  is  the  correct  statistical  factor  7  for  this  case, 

while  at  high  temperatures  it  leads  to  the  proper  saturation  value 

of  2N  -  N-. 

41  a 

The  solid  curves  in  Fig.  3  represent  the  fit  of  the  date  for 
the  three  lower  Te  pressures  using  Eq.  (3).  The  sample  fired  at 
the  highest  Te  pressure  waa  only  fitted  at  low  temperatures,  using 


i 
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Eq.  (1),  as  this  sample  was  not  purified  by  the  Zn  extraction 
method  and  was  known  to  contain  some  Cu.  In  Eq.  (3)  we  used 
the  previously  determined  value  of  0.048  eV  for  E^,  and  0.13  eV 
for  E 2,  although  a  range  of  values  (roughly  0.12  to  0.16  eV)  for 
E2  provided  reasonably  good  agreement.  The  acceptor  concentration*, 
Ng,  determined  from  the  fitting  of  the  data  for  the  samples  fired 
under  the  Te  pressures  of  10~^,  0.04  and  0.15  atm.,  were  approxi¬ 
mately  0.35  x  1017,  1.2  x  1017  and  1.7  x  10*7  cm  ^ ,  respectively. 
The  value  of  0.14  +  .02  eV  for  E2  is  admittedly  rather  close  to 
the  ionization  energy  for  Cu  acceptors.  The  high  temperature 
portion  of  the  curves  cannot,  however,  be  due  to  Cu  in  the  purified 
samples  for,  as  will  be  shown  below,  its  concentration  in  these 
samples  is  much  too  low.  It  is  of  course  possible  that  the  ZnTe 
crystals  contain  an  unidentified  acceptor,  unextrac table  by  the 
liquid  Zn  treatment,  in  concentrations  comparable  to  that  of  the 

Zn  vacancies,  and  that  the  observed  level  is  due  to  it.  We  believe 

21 

this  possibility  to  be  unlikely.  However,  before  the  deeper 
level  can  be  definitely  established  as  the  second  charged  state  of 
the  Zn  vacancy,  more  extensive  investigations  on  additional 
purified  and  analyzed  samples  over  a  more  extended  high  temperature 
range  would  be  required. 

Three  experiments  were  performed  to  assure  that  the  thermal 
treatment  to  which  the  crystals  were  subjected  during  the  diffusion 
process  was  not  producing  other  significant  changes  besides  the 
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incorporation  of  the  desired  impurities.  First,  Cu  was  diffused 

into  two  ZnTe  crystals  which  were  previously  fired  under  2  atm  Zn 

22 

pressure  and  0.5  atm  Te  pressure  respectively  and  which  exhibited 
the  0.048  eV  acceptor  levels.  In  each  case  the  incorporation  of  Cu 
resulted  in  the  disappearance  of  the  0.048  eV  level  and  the 
appearance  of  the  0.15  eV  level. 

In  a  complementary  experiment  two  crystals,  previously  doped 
with  Cu  and  exhibiting  the  0.15  eV  level,  were  fired  in  liquid  Zn 
and  liquid  Te,  respectively.  In  each  case  the  crystals  lost  their 
0.15  eV  levels  and  acquired  the  0.048  eV  levels.  One  crystal  was 
subjected  to  several  Cu  diffusion-Zn  extraction  cycles  and  was 
found  always  to  come  back  to  the  0.15  eV  level  after  Cu  diffusion 
and  to  the  0.048  eV  level  after  the  Zn  extraction.  An  analogous 
behavior  was  observed  with  an  Ag  diffusion-Zn  extraction  cycle 
and  an  Au  diffusion-Zn  extraction  cycle. 

64 

In  the  third  experiment  the  radioactive  isotopes  Cu  and 

Ag^°  were  separately  diffused  into  two  undoped  ZnTe  crystals. 

The  incorporation  of  the  radioactive  dopants  was  carried  out  by 

exactly  the  same  technique  as  that  used  with  the  crystal  on  which 

the  Hall  measurements  were  performed.  The  Cu  and  Ag  concentration 

after  the  diffusion,  as  determined  by  an  appropriate  counting 
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technique,  was  approximately  10  cm  .  By  several  successive 
etchings  and  countings  it  was  found  that  the  dopants  were  uniformly 
distributed  throughout  the  body  of  the  crystal.  After  extracting 
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both  crystals  with  liquid  Zn,  the  residual  concentration  of  Cu 
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and  Ag  in  them  was  found  to  be  approximately  10  cm  . 

These  experiments  show  that  the  introduction  of  the  employed 
acceptor  impurities  is  independent  of  whether  the  crystal  was 
previously  fired  under  Zn  or  Te  atmosphere.  Because  of  the 
reversibility  of  the  impurity  diffusion-Zn  extraction  cycles,  it 
is  also  certain  that  electrically  active  uncontrolled  impurities 
were  not  entering  the  crystals  in  significant  amounts  during  the 
various  firing  techniques  to  which  the  crystals  were  subjected. 

It  is  significant  to  note  that  the  analysis  of  the  carrier 
concentration  of  the  Cu  doped  sample  discussed  above  leads  to  an 
acceptor  concentration  in  fair  agreement  with  the  Cu  concentration 
determined  by  chemical  analysis.  The  experiment  with  the  radio¬ 
active  Cu  and  Ag  demonstrates  that  the  employed  diffusion  technique 
leads  to  a  uniform  distribution  of  the  dopant  in  the  crystal  in 
concentrations  where  it  exists  in  the  crystal  as  the  dominant 
impurity.  The  complementary  study  of  the  removal  of  the  Cu  and 
Ag  radioisotopes  by  Zn  extraction  shows  that  the  residual  con¬ 
centration  of  these  impurities  in  the  crystals  subjected  to  this 
purification  technique  is  several  orders  of  magnitude  below  the 
concentration  of  any  of  the  intentionally  introduced  acceptors 
which  have  determined  the  electrical  properties  of  the  crystals. 

Thus  we  see  there  are  acceptor  states  in  ZnTe  having  ionization 
energies  of  0.15*  0.11,  and  0.22  eV  definitely  associated  with  the 
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presence  of  Cu,  Ag,  and  Au  respectively.  The  simplest  inter¬ 
pretation  of  these  results  is  that  these  metal  impurities  go 
into  ZnTe  substitutionally  replacing  the  Zn.  The  acceptor  states 
are  then  directly  attributable  to  the  impurities.  One  feature  of 
this  picture  which  is  perhaps  disturbing  is  the  closeness  of  E 

0 

for  Cu  and  for  the  second  level  of  the  Zn  vacancy.  This  circumstance 
might  suggest  the  possibility  of  more  complex  models  for  the 
acceptors.  For  example,  it  might  be  assumed  that  Cu  acts  as  a 
donor  but  that  it  is  paried  with  the  vacancy.  The  0.048  eV  vacancy 
level  is  thereby  compensated.  While  complex  models  like  this 
cannot  be  completely  excluded  at  present,  their  complexity  alone 
makes  them  less  attractive.  Unless  further  evidence  is  adduced 
to  indicate  that  complexes  are  involved,  we  tend  to  favor  the 
simpler  model  of  the  acceptor.  The  near  equality  of  the  two 
ionisation  energies  is  believed  to  be  fortuitous. 

2.  ZnSe  -  Donor  Levels 

The  information  to  be  reported  below  about  donor  levels  in 
ZnSe  is  such  less  detailed  than  the  above  study  of  acceptor  levels 
in  ZnTe.  Although  the  Zn  extraction  method  was  found  to  remove 
acceptor  impurities  from  ZnSe  as  effectively  as  from  ZnTe,  as  yet 
little  is  known  of  its  effect  on  donor  impurities.  Therefore,  in 
the  experiments  with  n-type  ZnSe  to  be  described  below,  it  was  not 
possible  to  identify  the  donor  levels  with  particular  impurities 
or  native  defects. 
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Figure  4  a hows  the  results  of  Hall  measurements  for  a  ZnSe 
crystal  grown  with  lO1^  cm~^  Cl  (as  determined  by  chemical  analysis) 
and  subjected  to  one  extraction  with  liquid  Zn.  From  the  behavior 
of  between  100°K  and  room  temperature,  a  donor  level  at  0.19  eW 
is  determined.  This  value  is  very  close  to  the  0.21  eV  level 
found  by  Bube  and  Lind^  in  Br-doped  ZnSe.  It  is  possible  that 
these  levels  are  due  to  the  halogen  dopants  used,  although  the 
identification  has  yet  to  be  conclusively  established. 

The  previously  reported  donor  ionization  energies  in  ZnS  and 
ZnSe  are  in  the  range  of  0.2  to  0.3  eV  while  those  for  CdTe,  CdSe, 
and  CdS  are  much  smaller^  and  are  approximately  equal  to  those 
expected  on  the  basis  of  the  hydrogenic  model  (i.e.  0.01  to 


0.03  eV).  The  donor  depth  of  ZnSe  predicted  by  the  hydrogenic 

model,  Ed  *  13*6  (m*/m)  «~2  eV,  for  the  estimated  m*/m  of  0.15 

and  for  the  static  dielectric  constant  c  of  8.1,  is  0.03  eV. 

8 

This  is  almost  a  factor  of  ten  less  than  the  value  of  about 


0.2  eV  noted  above. 

The  question  arises  as  to  why  a  shallow  level  is  not  found 
in  ZnSe  (and  also  in  the  other  zinc  compounds)  while  shallow 
levels  are  found  in  the  cadmium  compound  semiconductors.  A 
simple  explanation  of  this  could  be  obtained  from  the  fact  that 
the  previous  observations  have  been  made  in  insulating  or  weakly 
n-type  crystals.  This  suggests  that  the  shallow  levels  were 
cos|>ietsly  compensated,  and  therefore  did  not  contribute  to  the 
free  carrier  concentration. 
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This  explanation  was  tested  by  preparing  a  more  strongly 

23 

n-type  sample  by  subjecting  an  undoped  crystal  to  two  successive 

Zn  extractions  in  an  attempt  to  reduce  the  concentration  of  the 

compensating  acceptors.  From  the  data  for  this  sample,  shown  as 

the  intermediate  curve  in  Fig.  4,  it  can  be  seen  that  the  carrier 

concentration  is  higher  and  the  slope  much  lower  than  for  the 

previous  sample.  Using  Eq.  (1)  with  g  *  2,  which  is  appropriate 

to  the  case  of  a  simple  donor,  the  data  are  fitted  using  - 

9.5  x  1016  cm”3,  Na  -  8.7  x  1016  cm”3  and  Ed  =  0.008  eV.  The 

calculated  points  are  shown  as  triangles  in  Fig.  4.  This  level  is 

indeed  shallow.  In  fact,  it  is  considerably  shallower  than  the 

hydrogenic  value.  At  these  moderately  large  values  of  and  Nft, 

the  latter  fact  undoubtedly  reflects  the  dependence  of  the  ionization 

energy  on  impurity  concentration  as  has  been  observed,  for  example, 
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in  Ge.  The  relatively  large  degree  of  compensation,  found  to 
be  present  in  this  sample  in  spite  of  the  extraction  of  acceptor 
impurities  and  the  suppression  of  the  concentration  of  native 
acceptors  (Zn  vacancies)  by  the  liquid  Zn  treatment,  demonstrates 
the  strong  tendency  of  ZnSe  toward  self -compensation. 

Finally,  a  sample  with  a  still  higher  carrier  concentration 

20  -3 

was  prepared  by  diffusing  1.5  x  10  cm  A1  into  a  previously 
purified  ZnSe  crystal.  The  Hall  constant  for  this  crystal  is 
given  as  the  lowest  curve  in  Fig.  4.  it  is  evident  from  the  fact 
that  Rjg  is  constant  from  about  16°K  to  370°K  that  the  ionisation 
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energy  has  vanished.  Evidently  the  inpurity  concentration  is 
greeter  then  the  critical  value  at  which  Ed  =  0,  and  the  Fermi 
level  either  coincides  with  the  conduction  band  edge  or  lies 
within  the  conduction  band.  The  sample  thus  represents  the  first 
case  of  degenerate  conduction  in  the  sine  II-VI  compounds. 

3.  Mobilities  of  ZnTe  and  ZnSe 

The  mobility  data  for  ZnSe  and  ZnTe  are  presented  in  Figs. 

3  and  6  respectively.  While  more  extensive  data  might  be  desirable 
for  a  detailed  analysis  of  the  mobilities,  the  present  results  are 
quite  adequate  to  permit  a  study  of  the  dominant  scattering 
mechanisms.  From  the  near  equality  of  the  mobilities  (u)  for 
T  >  100°K  for  samples  with  significantly  different  defect  concen¬ 
trations  and  from  their  fairly  rapid  increase  with  decreasing 
temperature,  it  appears  that  the  ix  for  T  >  100°K  are  determined 
by  the  intrinsic  properties  of  the  two  materials  and  not  by  the 
crystal  defects.  In  addition  it  will  be  shown  below  that  the 
magnitude  and  temperature  dependence  can  in  fact  be  accounted  for 
by  a  lattice  scattering  mechanism.  Below  about  100°K  on  the  other 
hand,  charged  impurity  scattering  is  indicated  by  the  increase  in 
U  with  temperature.  This  is  supported  by  the  study  of  u  for  the 
one  sample  for  which  low  temperature  data  is  available  (lower 
experimental  curve  in  Fig.  3) •  Using  the  donor  and  compensating 
acceptor  concentrations  obtained  from  the  analysis  of  the  carrier 
concentration  and  the  Brooks -Herring  formula,  we  obtain  the  dotted 
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curve  in  Fig.  5»  The  agreement  Is  satisfactory  considering  the 
uncertainties  in  the  data  and  in  the  single  application  of  the 
theory. 

It  should  be  noted  that  a  study  of  the  mobilities  of  ZnTe 
and  ZnSe  will  in  a  general  way  parallel  the  corresponding  studies 

of  the  more  thoroughly  investigated  semiconductors  in  the  III-V 

25 

compound  family.  For  as  Ehrenreich  has  shown,  at  least  for  the 
n-type  materials,  the  polar  character  of  these  semiconductors 
plays  a  major  role  in  their  transport  properties.  This  would 
also  be  expected  to  be  true  for  the  II-VI  compounds  because  of 
their  strongly  polar  nature.  To  provide  a  reasonably  complete 
study  of  the  intrinsic  (or  lattice)  mobility  of  these  materials, 
four  scattering  mechanisms  should  be  considered.  These  are 
scattering  by  the  piezoelectric  activity  of  the  acoustic  modes, 
the  deformation  potential,  and  the  non -polar  and  polar  Interactions 
with  the  optical  modes.  For  all  four,  u  drops  off  with  increasing 
T.  As  we  will  attempt  to  determine  which  of  these  mechanisms 
contribute  apprecialby  to  the  observed  u,  we  will  in  the  following 
consider  each  of  them  separately  and  in  some  detail.  Unless 
otherwise  indicated,  we  will  avoid  the  complications  arising  from 
the  degenerate  valence  bands  (for  p-type  ZnTe)  by  the  relatively 
rough  approximation  of  treating  them  as  a  simple  band.  With  a 
suitable  m*  this  should  provide  us  with  results  which  are  at 
least  semi -quantitatively  correct;  and  this  is  all  that  will  be 
required  for  our  purposes. 
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The  scattering  by  piezoelectrically  active  acoustic  nodes, 

2&  27 

first  treated  by  Meijer  and  Polder  and  Harrison,  results  fron 
the  snail  component  of  polarization  present  in  the  long  wavelength 
acoustic  nodes.  The  nobility  of  carriers  in  a  sluple  band  with, 
this  interaction  is 

upiezo  "  1,05  p  <ul>€s  el4  T~1/2  c“2/v  8ec*  (4) 

where  eg  is  the  static  dielectric  constant,  e^  the  piezoelectric 
constant  (in  esu/cn  ),  p  the  density,  and  <u^>  the  square  of 
longitudinal  sound  velocity  averaged  over  direction.  Aside  fron 
n*/m,  for  which  we  use  the  values  stated  above,  the  values  of  the 
various  parameters  used  in  Eq.  (4)  are  given  in  Table  1.  The 
resulting  nobilities  are: 

Upiezo^ZnTe^1’6  x  1q5  <  300/n1/2  cn2/V  sec. 

Upie20(znse)«3.5  x  105  (300/T)1/2  cm2/V  sec. 

These  values  exceed  the  measured  ones  by  well  over  a  factor  of  100 
and  thus  Indicate  that  this  nechanisn  does  not  play  an  important 
role  in  determining  the  u  for  either  of  the  two  materials. 

The  theory  of  deformation  potential  (or  acoustic  mode) 
scattering  is  wall  known.  Tha  mobility  for  this  casa  can  be 
expressed  as 

udp  "  3,0  x  l0"5  (**/*) 5/2  P  <“j>  t*3/2  Et>2  c®2/v  MC  (5) 

where  denotes  tha  deformation  potential  (in  eV)  for  tha  relevant 

band  edge.2® 
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Unfortunately  there  Is  no  direct  information  about  these 
quantities  for  ZnSe  and  ZnTe.  However,  information  about  the 
pressure  coefficients  of  a  fairly  large  number  of  semiconductors 
is  available.  From  the  regularity  of  the  pertinent  data  it  appears 
that  estimates  of  the  required  quantities  can  be  made  which  are 
sufficiently  accurate  for  our  purpose.  We  take  the  values  of  4  ev 
and  2  ev  as  the  deformation  potentials  for  the  conduction  band  edge 

2Q 

of  ZnSe  and  the  valence  band  edge  of  ZnTe  respectively.  With 
these  values  the  calculated  mobilities  are 

Udp(ZnTe)  -  4  x  103  (300/T)1/2  cm2/V  sec 

4  1  /o  o 

Udp(ZnTe)  -  2  x  10  (300/T)  '  cm  /V  sec. 

Since  these  values  are  about  a  factor  of  40  times  larger  than  the 
measured  values,  it  appears  that  acoustic  mode  scattering  is  not 
important  for  these  compounds.  Of  course,  if  the  deformation 
potentials  are  much  larger  than  our  estimated  values,  (say  by  a 
factor  of  five)  our  conclusion  would  be  in  error.  However,  it 
seems  rather  unlikely  in  light  of  present  knowledge  that  these 
quantities  can  be  so  huge. 

The  expression  for  the  mobility  of  holes  in  degenerate  p-llke 

bands  scattered  by  the  non-polar  Interaction  with  the  optical 

phonons,  u. _ ,  can  be  obtained  from  the  work  of  Ehrenreich  and 

npo 

Overhauser.30  From  their  study  of  the  mobility  of  holee  in  Ge, 


it  c#n  be  readily  shown  that 


1.35  X  1017  PQ^g3  r 

.57?  075 
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x  •»p(-8T~1x)d« 
~  ?(*"1) 


£/f 


with  f  (t)  -  n(l+t)1/2  +  (n+l)(l-t)1/2  for  t  <  1 


n(l+t) 


1/2 


for  t  >  1. 


(6) 


Here  n  Is  the  usual  optical  phonon  occupation  number  given  by 

[exp(0/T)  -  l]-1,  9  the  Debye  temperature,  m,  and  m2  the  meases  of 

the  light  and  heavy  hole  respectively,  a  the  lattice  constant  (in 

cm),  and  the  coupling  constant  (in  eV).  The  constant  is 

analogous  (but  not  equal)  to  the  deformation  potential  for  the 

31 

optical  modes  in  say  Conwell's  formulation.  In  line  with  our 
limited  knowledge  about  the  masses  for  ZnTe  the  factor  involving 

*“5/2  _c/p 

the  masses  is  replaced  by  m  ■  (0.6m)  .  The  biggest 

uncertainty,  however,  is  in  the  magnitude  of  the  optical  mode 
coupling  parameter  C^,  about  which  little  is  known.  To  our 
knowledge  the  only  semiconductor  for  which  the  strength  of  this 

coupling  has  been  determined  is  Ge.  In  their  recent  study  of 

32 

the  mobility  of  p-type  Ge,  Brown  and  Bray  have  determined 
coupling  parameters  for  the  different  modes  (so  called  optical 
and  acoustic  deformation  potentials).  From  their  results  it  is 
possible  to  deduce  that  the  value  of  is  20  eV,  which  corresponds 
to  their  optical  mode  deformation  potential  of  8.8  eV.  On  the 
other  hend  it  is  possible  to  meke  a  rough  estimate  of  the  upper 


erne 
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limit  of  this  parameter  for  InSb.  We  find  the  value  to  be  15  eV.^ 

Using  the  value  of  15  eV  for  ZnTe  we  find  that  u.  has  the  values 

npo 

of  790,  2000,  and  1800  cm2/V  sec  at  T  -  300°,  200°,  and  100°K 
respectively.  These  values  arm  too  large  by  a  factor  of  about 
8  at  300°K  and  25  at  100°K. 


With  significantly  larger  values  for  the  mass  and  than  the 
estimates  that  we  have  used,  ^npQ  could  yield  values  comparable 
to  the  observed  values.  On  the  other  hand  it  is  important  to  note 
that  the  temperature  dependence  of  M-np0  also  is  not  in  very  good 
accord  with  the  data,  as  it  increases  too  rapidly  with  decreasing 
temperature.  Thus  while  we  cannot  at  present  definitely  rule 
out  the  possibility  that  this  type  of  scattering  contributes 
significantly  to  the  total,  these  observations  suggest  that  it 
doesn't  dominate  the  scattering. 

This  mechanism  can  be  completely  ruled  out  for  n-type  ZnSe 
for  two  reasons.  First,  the  effective  mass  is  small  and  u  oC  (m*) 
The  second  and  more  important  reason  is  that  for  the  (s-like) 
conduction  band  minimum  the  matrix  element  for  scattering  between 
electronic  states  k  and  k'  vanishes  to  lowest  order  in  the  phonon 
wave  vector  q  =  k  -  k ' .  The  reduction  in  the  scattering  due  to 
this  selection  rule  is  ~  10 

The  last  mechanism  to  be  considered  is  the  scattering  of  the 
carriers  by  the  electric  polarisation  associated  with  the  optical 


-5/2 
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■odes.  The  strength  of  the  interaction  is  indicated  by  the  polar 
coupling  constant  a  which  is 

a  -  (tn^/m)1^2  (Ry/fla^)1^2  (e^1  -  e"1) 


where  Ry  is  the  Rydberg  of  energy,  the  energy  of  the  longitudinal 
optical  mode  for  long  wavelength,  and  c  the  high  frequency  or 

-r 
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optical  dielectric  constant.  Using  the  relevant  parameters,  which 

£ 

are  given  in  Table  I  and  the  assumed  values  of  m  / m  we  find  that 
a  -  0.41  and  0.39  for  electrons  in  ZnSe  and  holes  in  ZnTe  respectively. 
By  noting  that  the  factor  which  more  properly  reflects  the  magnitude 
of  succeeding  orders  of  perturbation  is  a/6,  it  is  evident  that 
the  weak  coupling  approach  is  probably  adequate  for  these  materials. 
The  perturbation  theory  result  for  simple  parabolic  bands  obtained 
by  Howarth  and  Sondheimer^  and  in  which  the  Callen  effective  charge^ 
is  used  can  be  written  as 


upolar 


0.870(a/tn*i 


]  G(z)e 


-5 


cm  /V  sec 


(7) 


where  z  =  fku^/kT  =  d/T,  is  in  eV,  and  G(z)e  is  a  tabulated 

function.  A  more  general  result  having  the  same  form  as  Eq.  (7) 

but  in  which  the  screening  by  the  carriers  is  taken  into  account 

37 

was  given  by  Ehrenreich.  His  G(z)  is  a  function  of  the  carrier 

2 

concentration  through  the  plasma  frequency  co  defined  by  u>  * 

P  P 

4ve2n/«*.  The  difference  between  the  G(z)  for  the  screened  and 
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unscreened  interactions  vanishes  as  to  /cu,  — >  0.  Since  to  /o>,  «  1 

P  •  P  * 

for  all  ZnSe  and  ZnTe  crystals  considered  in  this  context,  the 
screening  can  be  neglected. 

The  values  obtained  from  Eq.  (7)  are  shown  as  the  solid  curves 
in  Figs.  5  and  6.  in  contrast  to  the  previous  mechanisms  considered, 
polar  scattering  leads  to  results  which  are  quite  close  to  the 
experimental  mobilities.  It  is  noteworthy  that  there  are  no 
adjustable  parameters  in  M-^Q^ar  except  perhaps  the  m*s  (which,  in 
fact,  have  not  been  varied).  There  is  an  apparent  small  discrepancy 
in  that  the  calculated  curves  exhibit  a  small  positive  curvature 
at  the  higher  temperatures  while  the  data  do  not.  However,  as 
has  been  found  in  a  recent  study  of  CdTe"5  this  can  be  simply 
understood  in  terms  of  the  temperature  dependence  of  the  static 
dielectric  constant  which  has  been  neglected  in  the  solid  curves. 

A  small  variation  in  c  (T)  has  an  appreciable  effect  on  a  because 

9 

of  the  considerable  cancellation  between  e  1  and  e  To  estimate 

co  s 

this  effect  we  have  made  the  assumption  that  the  temperature 

variation  of  e  for  ZnTe  and  ZnSe  is  similar  to  that  for  CdTe. 

8 

It  is  to  be  noted  that  the  dashed  curves  which  include  this 
correction,  have  a  temperature  dependence  which  more  closely 
resembles  the  data  in  the  intrinsic  range. 

In  summary,  the  polar  optical  mode  scattering  clearly 
dominates  the  intrinsic  scattering  of  electrons  in  ZnSe,  the 
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scattering  frost  all  other  mechanism  being  very  nuch  weaker.  This 
mechanism  is  undoubtedly  also  very  important  for  p-type  ZnTe  as 
the  estimated  mobility  due  to  it  alone  has  the  correct  magnitude 
and  temperature  dependence.  The  other  scattering  mechanism  that 
possibly  could  be  important  in  the  scattering  of  the  holes  is  the 
non-polar  interaction  with  the  optical  phonons.  This  could  be  the 
case  if  the  appropriate  coupling  parameter  is  considerably  larger 
than  the  values  we  have  considered.  Finally,  it  should  be  noted 

that  high  temperature  data  would  be  very  helpful  in  assessing  the 

o/o 

roles  of  these  two  mechanisms  since  u  oc(9/T)  while 

npo 

^polar 06  (0/T)1/2  £or  1  »  0- 

IV.  CONCLUSION 

The  present  study  of  the  electrical  transport  properties  of 
ZnSe  and  ZnTe  has  shown  that  shallow  impurity  states,  observed 
in  materials  like  Si,  Ge  and  the  more  thoroughly  studied  compounds 
of  the  III-V  compound  family  are  more  widespread  in  the  II-VI 
compounds  than  has  been  believed  in  the  past.  Furthermore,  it 
was  found  that  the  processes  controlling  the  electircal  conduction 
of  the  investigated  compounds  are  similar  to  those  that  have  been 
shown  to  be  Important  in  such  compounds  as  InSb  and  GaAs.  Thus, 
in  spite  of  their  less  covalent  character,  the  II-VI  compounds 
are  in  these  respects  quite  similar  to  the  more  conventional 
sesdconductors. 
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The  Method  of  firing  the  crystals  in  liquid  Zn  has  enabled 
us,  by  sufficiently  suppressing  the  concentration  of  compensating 
impurity  and  native  defects,  to  observe  a  shallow  donor  level  in 
ZnSe  with  ionization  energy  of  ~  0.01  eV.  It  has  also  proved 
possible  to  prepare  heavily  Al-doped  ZnSe  crystals  which  are 
degenerate . 

The  discovery  of  shallow  acceptors  in  ZnTe  represents  the 
first  case  in  which  uncompensated  acceptor  levels  of  near  hydro- 
genic  depths  have  definitely  been  observed  in  any  of  the  II-VI 
semiconductors.  Acceptor  levels  correlated  with  Cu,  Ag,  and  Au 
dopants  were  identified  and  found  to  have  ionization  energies 
of  0.15,  0.11  and  0.22  eV  respectively.  The  levels  are  most 
probably  simple  acceptors  resulting  from  the  substitutional 
replacement  of  Zn  atoms  by  noble  metal  atoms;  although  more 
complex  models  cannot  presently  be  excluded.  A  level  at  0.048 
eV  in  undoped  crystals  has  been  identified  as  being  due  to  a 
Zn  vacancy.  Some  evidence  is  presented  for  what  is  believed  to 
be  the  second  charge  state  of  the  vacancy. 

The  Hall  mobility  data  on  lightly  doped  ZnSe  and  ZnTe  strongly 
suggest  that  the  mobilities  for  T  >  100°K  are  determined  by  the 
intrinsic  properties  of  the  crystals.  Analysis  of  the  intrinsic 
(or  lattice)  mobility  indicates  that  in  the  case  of  ZnSe  polar 
optical  mode  scattering  strongly  dominates  over  the  scattering 
by  the  piezoelectric  activity  of  the  acoustic  modes,  the  deformation 
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potential,  and  the  non-polar  interaction  with  the  optical  nodes. 
Similarly,  the  polar  scattering  of  the  optical  phonons  appears 
to  be  the  inportant  mechanism  determining  the  mobility  of  holes 
in  ZnTe.  It  is  possible,  however,  that  the  non-polar  optical 
mode  interaction  could  also  contribute  appreciably  to  the  scattering 
of  the  holes  if  the  relevant  coupling  parameter  is  significantly 
larger  than  the  corresponding  parameters  for  Ge  end  InSb. 
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TABLE  I 

Parameters  Used  in  Calculating 
Mobilities  for  ZnSe  and  ZnTa 


ZnSe 

ZnTe 

«#(T  -  300°K) 

8.1* 

10. 1C 

e 

00 

5.75a 

8.26c 

4ia>j(eV) 

0.0314* 

0.0259b 

/Stat.  coul • x 

*14'  2  > 

cm 

1.5  x  104  d 

8.5  x  lo3  d 

p  «*i>  (*¥) 

1.06  x  1012  * 

8.5  x  10U  * 

cm 


a)  reference  4 

b)  reference  11 

c)  reference  1 

d)  D.  Berl incourt,  H.  Jaffe  and  L.  1.  Shioaawa,  submitted 
for  publication  in  Phys.  Rev. 

a)  derived  from  the  elastic  constants  in  d. 
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FIGURE  CAPTIONS 

Fig.  1  Temperature  dependence  of  the  Hall  coefficient  for  ZnTe 
crystal#  doped  with  Au  (squares),  Cu  (triangles)  and  Ag 
(circles)  and  for  an  undoped  ZnTe  crystal  fired  in  liquid 
Zn  at  900°C  (diamonds). 

Fig.  2  Temperature  dependence  of  the  free  hole  concentration  for 

a  ZnTe  crystal  found  (by  chemical  analysis)  to  contain 
16 

(2  +  1)  x  10  cm"*5  Cu.  The  data  are  indicated  by  circles. 

The  solid  line  represents  the  hole  concentration  calculated 
from  Eq.  (1)  with  Nft  -  3-5  *  1016  cm'3,  Nd  -  3*1  *  1015  cm'3 
and  E  -  0.149  eV. 

O 

Fig.  3  Temperature  dependence  of  the  hole  concentrations  for 

undoped  ZnTe  crystals  fired  under  various  partial  pressures 
of  Te.  Theoretical  curves,  calculated  from  Eq.  (3),  are 
indicated  by  solid  lines.  The  dashed  curve  represents 
the  low  temperature  range  for  the  pT#  -0.26  atm  crystal 
fitted  using  Eq.  (1)  (see  text). 

Fig.  4  Temperature  dependence  of  the  Hall  coefficient  for  ZnSe 
crystals  doped  with  Cl  (diamonds)  and  A1  (squares)  and 
for  an  undoped  ZnSe  crystal  (circles).  The  points 
calculated  from  Eq.  (1)  for  the  undoped  crystal,  with 
Nd  -  9*5  x  1016  cm*3,  N0  ■  8.7  x  lO1^  cm'3  and  Ed  -  0.00ft  e?, 
are  indicated  by  triangles. 
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fit*  5  Temperature  dependence  of  the  Hall  mobility  for  an  undopad 

16  “3 

ZnSe  crystal  with  n30o°K  ■  6.9  x  10  cm  3  (squares),  a 
Cl-doped  crystal  with  n30o°K  “  x  c®  ^  (triangles) 

and  a  Cl-doped  crystal  with  n3Q0°£  -  6.8  x  10  cm  3 
(circles).  The  solid  curve  represents  the  calculated 
mobility  for  polar  scattering  of  optical  phonons.  The 
dashed  curve  includes  an  approximate  correction  for  the 
temperature  variation  of  €a«  The  dotted  curve  is  the 
mobility  calculated  for  charged  impurity  scattering  for 
the  undoped  sample. 

1 

Fig*  6  Temperature  dependence  of  the  Hall  mobility  for  two 

undoped  (circles  and  squares)  and  a  Ag-doped  (triangles) 

ZnTe  crystal.  The  solid  curve  represents  the  calculated 

1 

mobility  for  polar  scattering  of  optical  phonons.  The 
dashed  curve  includes  an  approximate  correction  for  the 
temperature  variation  of  «a • 
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C.  Light  Emission  From  GaAs^_xPx  and  GaSb  Junctions  (J.  H.  Racetta) 

Work  Is  continuing  on  the  study  of  light -emit ting  GaAs^_xPx 
p-n  junctions  using  crystals  prepared  by  growth  from  a  gallium 
solution,  and  by  halogen  transport.  The  latter  crystals  were 
made  available  to  us  from  another  source  within  the  laboratory. 

It  is  hoped  that  these  studies  will  help  to  improve  the  reproducibility 
and  luminescence  efficiency  of  these  mixed  crystal  junctions,  which 
are  at  present  unsatisfactory. 

Diodes  with  very  good  electrical  characteristics  have  been 
made  from  zinc  diffused  (  96o°C,  1  hour)  n-type  material. 

Structures  with  carefully  polished  sides  have  been  made  with  the 
hope  of  obtaining  laser  action,  with  consequent  coherence  of  the 
visible  junction  radiation.  As  yet,  however,  sufficiently  planar 
junctions  have  not  been  obtained.  Presumably  this  is  due  to  the 
poor  crystal  perfection  of  the  n-type  bulk  material  available  thus 
far.  The  exact  location  and  shape  of  the  junction  have  been 
determined  either  by  viewing  the  radiation  produced  at  the  junction, 
or  by  preferential  etching  of  the  p-type  material  electrolytically 
in  a  25°/o  NaOH  solution.  Capacitance  measurements  indicate  that 
the  junctions,  apart  from  non-planarity,  are  somewhat  thinner 
than  similarly  fabricated  GaAs  junctions,  being  typically  300X 
wide  at  zero  bias. 

The  light  output  of  diodes  fabricated  here  has  been  found  to 
vary  linear ly  with  diode  current.  The  spectral  distribution  of 
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this  radiation  has  been  measured  as  a  function  of  current  at  77°K 
by  J.  Kingsley.  An  appreciable  shift  to  lower  wavelength  of  the 

peak  from  7500X  to  6800X  was  observed  as  the  current  was  increased 

2 

from  8  to  800  ma  for  a  junction  of  0.005  cm  area.  No  line 
narrowing  was  noted,  the  width  at  half  maximum  of  the  energy  peak 
being  around  3008. 

Work  is  also  being  carried  out  on  GaSb  as  a  radiation  source. 
This  is  another  "direct"  material  which  consequently  should  yield 
high  efficiency  junction  recomb ination  radiation. 

The  procedure  used  here  is  the  same  as  that  used  in  connection 
with  GaAs^P^^.  Diodes  with  polished  sides  have  been  fabricated 
from  n-type  GaSb  into  which  zinc  has  been  diffused  at  66o°C  for 
90  hours. 

The  radiation  from  such  a  junction  at  77°K  has  been  measured 
by  W.  E.  Engeler.  A  shift  of  the  radiation  peak  from  1.75u>  to  1.70u 

was  observed  as  the  diode  current  was  varied  from  15  ma  to  150  ma, 

2 

the  junction  area  being  0.004  cm  .  Further  increase  in  the  current 
caused  only  a  very  slight  decrease  in  wavelength  of  the  peak.  The 
width  of  this  peak  at  half  maximum  was  found  to  be  about  0.17u> 

The  radiative  output  varied  linearly  with  current  over  most  of 
the  range  explored,  though  it  increased  at  a  lower  rate  for  low 
currents.  More  extensive  investigation  of  this  material  is  planned. 
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